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PART  A 


The  work  of  G.  F.  Carrier  under  these  contracts  can  best  be  sep- 
arated in  two  parts: 

Part  I (research  which  has  been  published  in  the  papers  listed  below) 

Under  Contract  N00014-67-A-0298-0002 

"Singular  Perturbation  Theory  and  Geophysics,"  Von  Neumann  lecture  given 
at  the  American  Mathematical  Society  Meeting,  Aug.  23-29,  1969. 

SIAM  Review,  Vol.  12,  No.  2,  April  1970. 

"Stochastically  Driven  Dynamical  Systems,"  J.  Fluid  Mech.,  44,  part  2,  pp. 
249-264,  1970. 

"A  Model  of  a Mature  Hurricane,"  (with  A.  L.  Hammond  and  0.  D.  George), 

J.  Fluid  Mech..  Vol.  47,  part  1,  pp.  145-170,  1971. 

"The  Dynamics  of  Tsunamis,"  American  Math.  Soc.  1970.  Also  In  Mathematical 
Problems  in  the  Geophysical  Sciences.  Vol.  1 - Geophysical  Fluid 
Dyrtainics . William  H.  Reid,  Editor,  1971. 

"Another  Useful  Singular  Perturbation  Device,"  Congres  International  des 
Mathematlciens,  Nice,  Sept.  1970. 

Under  Contract  N00014-67-A-0298-0033 

"The  Effects  of  Strain  Rate  on  Diffusion  Flames,"  (with  F.  E.  Fendell  and 
F.  Marble),  SIAM  J.  AppI.  Math..  Vol.  28,  No.  2,  pp.  436-500, 

March  1975. 

"The  Effect  of  Strain  Rate  on  Diffusion  Flames.  Part  11.  Large  Straining," 
(with  F.  E.  Fendell),  SIAM  J.  Anol.  Math..  Vol.  30,  No.  3,  pp.  515-527, 
May  1976. 


□ □ 
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Part  II.  This  part  consists  of  several  studies  which  were  Initiated  under 
N00014-67-A-0298-0033  (some  completed),  which  have  not  yet  been  published: 

(a)  Wave  Attenuation  In  Harbors.  This  Is  an  effort  to  understand  the 
long  ringing  time  of  Tsunami  Instigated  waves  In  Inlets  and  harbors,  par- 
ticularly the  long  (3-day)  ringing  of  the  S-mlnute  slashing  mode  In  Nawlllwlll 
Bay  just  after  the  Alaskan  earthquake.  This  work  Is  being  continued  under 
other  sponsorship  by  Claude  Nolseaux,  a graduate  student,  with  the  supervision 
of  Carrier. 

(b)  Wave  Reflection  by  Irregular  Topography.  This  completed  effort 
provided  a reasonably  efficient  way  of  calculating  the  reflection  of  gravity 
waves  by  an  Irregular  bottom  contour.  It  also  arose  In  connection  with 
Tsunami  propagation  but,  until  a specific  application  seems  to  merit  a detailed 
application,  publication  Is  not  likely  to  arouse  much  Interest. 

(c)  The  Mid-Latitude  Ocean  Circulation.  For  many  years,  those  who 
study  models  of  the  wind-driven,  mid-latitude  circulation  have  failed  to  treat 
successfully  the  steady  (but  unstable)  state  corresponding  to  the  must  Inter- 
esting Rossby  number  range.  A new  formulation  and  an  appropriate  (classical) 
technique  have  yielded  very  encouraging  preliminary  results  and  this  work  Is 
continuing  under  other  sponsorship. 

(d)  During  a sabbatical  visit  to  Harvard  by  Wilbert  Lick  of  Case-Western 
Reserve,  he  and  Carrier  tried  to  model  the  dynamics,  thermo-dynamics  and  heat 
transfer  mechanisms  which  are  Important  in  Cumulo-Hlmbus  activity.  The  work 
was  Instructive  but  no  definitive  results  were  obtained.  The  Insight  gained 
has  been  useful  In  other  studies  but  no  publications  of  this  work.  In  Its  own 
right,  seemed  justified. 

(e)  During  the  spring  term  of  1971-72,  Leon  Knopoff 's  visit  to  Harvard 
had  ONR  support  (letter  report  by  G.  F.  Carrier  of  June  15,  1973).  His  pub- 


llcaClon  "The  Dynamics  of  a One-Dimensional  Fault  in  the  Presence  of 
Friction,"  Publication  No.  1166,  Institute  of  Geophysics  and  Planetary 
Physics,  University  of  California,  Los  Angeles  was  an  outgrowth  of  that  visit 


PART  B 


The  work  of  F.  H.  Abernathy  under  this  contract  can  be  divided  into 
separable  though  Interconnected  activities.  The  first  activity  was  the 
development  of  a technique  for  using  laser -Doppler  anemometry  in  high 
strain  rates  flows.  The  second  activity  was  a detailed  investigation  of 
the  mixing  in  a turbulent  boundary  using  drag-reducing  polymers  to  selectively 
inhibit  motions  based  on  their  strain  rates.  Both  of  these  activities  were 
initially  supported  under  this  contract,  with  continuing  sponsorship  provided 
by  NSF. 

The  work  performed  under  both  activities  has  been  published: 

(1)  "Turbulence  Spectra  Using  Laser-Doppler  Anemometry  and  Selective 

Seeding,"  F.  H.  Abernathy,  J.  R.  Bertschy,  and  R.  W.  Chin, 
Biennial  Turbulence  Conference,  University  of  Missourl- 
Rolla,  Oct.  1977. 

(2)  "Modifications  to  Laminar  and  Turbulent  Boundary  Layers  Due  to 

the  Addition  of  Dilute  Polymer  Solutions,"  J.  R.  Bertschy 
and  F.  H.  Abernathy,  BHRA  Fluid  Engineering,  Second  Inter- 
national Conference  on  Drag  Reduction,  Cambridge,  England, 
August  1977. 

Both  of  these  papers  will  be  available  In  the  proceedings  of  these 
meetings.  For  convenience,  a copy  of  each  paper  is  attached  to  this  report 
in  Appendix  I-a,  and  Appendix  1-b,  and  briefly  summarized  as  follows: 


6 


Using  the  selective  seeding  technique  described  In  (1)  above.  It  Is 
possible  to  obtain  thousands  of  velocity  determinations  per  seconds  In 
water  shear  flows  with  spatial  resolution  of  the  order  of  25  M.  A plane 
of  hydrogen  bubbles  Is  generated  by  electrolysis  from  a wire  upstream 
of  the  laser-Doppler  scattering  volume.  The  Intersection  of  the  plane 
of  bubbbles  with  the  scattering  volume  becomes  the  major  source  of 
scattered  light  If  the  flow  Is  prefiltered.  This  technique  requires  a 
relatively  simple  correction  for  velocity  determinations.  The  correction 
Is  necessary  because  the  bubbles  are  in  the  wake  of  the  generating  wire. 

The  bubble  generation  rate  can  be  controlled.  At  the  upper  limit  a 
hydrogen  bubble  Is  nearly  always  crossing  the  scattering  volume.  The 
laser-Doppler  output  then  Is  nearly  a continuous  time  signal.  Hence  the 
laser-Doppler  signal  can  be  time  averaged,  frequency  filtered,  etc. 

Power  spectural  density  measurements  can  be  made  In  flows  where  It  would 
have  been  Impossible  before  the  technique  was  developed.  It  also  makes 
it  possible  to  obtain  detailed  velocity  profiles  of  flows  which  would 
have  been  Impossible  to  obtain  by  other  means.  It  Is  felt  that  this 
technique  will  become  useful  In  a variety  of  flow  situations. 

The  addition  of  polymers  In  dilute  concentrations  to  laminar  and 
turbulent  boundary  layer  causes  significant  effects.  What  has  been 
called  "early  turbulence"  by  others  has  been  shown  to  be  not  turbulence 
but  rather  fluctuations  In  the  laminar  layer,  probably  associated  with  the 
relaxation  time  of  the  polymer  Itself.  Using  the  selective  seeding 
technique,  it  has  been  possible  to  analyse  the  spectrum  of  the  fluctiutions 
and  show  that  they  are  distinct  from  turbulence.  In  addition  as  polymer 


concentration  and  strain  rate  Increase,  the  difference  between  laminar 
and  turbulent  layers  disappears.  Implications  of  these  results  are  under 
active  Investigation. 
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PART  C 

Under  this  contract  a substantial  research  effort  on  "Grid  Turbulence 
In  a Liquid  System"  was  undertaken  by  Carol  Russo  under  Professor 
Howard  W.  Emmons. 

The  final  work  on  this  experiment  was  completed  In  May  of  1975. 

At  the  same  time  a Ph.  D.  thesis  was  written  and  accepted  entitled, 

"An  Experimental  Study  of  Grid  Turbulence  in  a Liquid".  A copy  of  the 
thesis  synopsis  Is  attached  which  summarizes  the  entire  effort  on  this 
project  and  draws  conclusions. 


j 
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SYNOPSIS 

Measurements  were  made  of  the  Intensities  and  of  the  transverse  and 
longitudinal  double  correlations  of  grid  turbulence  In  the  Initial  region 
of  decay  in  methanol  to  verify  for  a liquid  the  theoretical  and  seml- 
emplrlcal  equations  for  isotropic  turbulence  that  have  been  previously 
verified  In  air.  This  experiment  was  undertaken  primarily  for  two  reasons. 
First,  some  question  of  the  validity  of  these  equations  for  a liquid  was 
raised  by  Tan  and  Ling's  Intensity  data  taken  behind  a grid  in  water  because 
their  Intensities  were  much  higher  than  those  measured  in  air  under  apparently 
equivalent  experimental  conditions.  It  was  not  clear  whether  the  differences 
in  molecular  structure  between  gases  and  liquids  which  change  the  mechanisms 
Involved  In  diffusion  processes  would  effect  the  basic  equations  describing 
the  behavior  of  isotropic  turbulence.  Secondly,  it  became  clear  in  the  air 
data  that  changes  In  the  mean  speed  U and  the  mesh  size  M which  partly 
comprise  the  mesh  Reynolds'  number  (UM/v),  the  primary  similarity  parameter, 
did  not  equally  effect  the  resulting  intensity.  It  seemed  important  to 
determine  If  changing  the  kinematic  viscosity  would  significantly  effect 
the  properties  of  the  grid  turbulence.  The  experiments  in  air  could  not 
vary  the  kinematic  viscosity  over  a large  range.  Using  a liquid  Instead  of 
a gas  Is  a simple  way  to  dramatically  vary  the  kinematic  viscosity. 

Measurements  of  both  Intensity  and  correlation  were  in  general  In  good 
agreement  with  the  experimental  results  found  in  air.  The  grid  turbulent 
energy  decayed  linearly  with  distance  with  slopes  equal  to  those  in  air. 

The  turbulent  intensity  reacted  to  changes  in  similarity  parameters  such  as 
■esh  to  wire  diameter  ratio  and  mesh  and  wire  diameter  Reynolds'  numbers 
In  a manner  similar  to  that  found  In  air.  Also,  the  Intensities  In  methanol 
fell  within  the  zone  of  scatter  of  the  air  data.  However,  detailed  ccoparlsons 


of  the  values  of  the  grid  turbulence  intensities  show  differences  much 
larger  than  known  experimental  errors  and  which  cannot  be  explained  by 
the  differences  in  factors  such  as  the  free  stream  turbulence.  For  the 
same  values  of  the  similarity  parameters  Re^,  M/d,  and  x/M,  the  intensities 
in  methanol  were  lower  than  those  found  in  air.  The  lateral  and  longitudinal 
double  correlations  were  in  good  agreement  with  those  measured  in  air  also. 
The  correlations  measured  in  methanol  were  also  in  agreement  with  the  purely 
theoretical  equation  relating  the  lateral  and  longitudinal  correlations  for 
Isotropic  turbulence  which  shows  that  in  methanol  grid  turbulence  is  also 
Isotropic.  In  contrast  to  the  grid  intensities,  the  behavior  of  the 
correlations  is  well  explained  in  detail  by  the  usual  similarity  parameters 
in  both  methanol  and  air. 

It  is  clear  that  while  the  general  behavior  of  grid  turbulence  has  been 
given  by  the  theoretical  and  seml-emplrlcal  equations  verified  in  air,  some 
of  the  parameters  needed  to  predict  the  value  of  the  intensity  precisely 
remain  to  be  found.  It  would  require  an  experiment  much  more  extensive  and 
with  smaller  experimental  errors  to  unravel  the  details  of  grid  turbulence 


behavior . 
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MCX)IFICATION5  TO  LAMINAR  AND  TURBULENT  BOUNDARY  LAYERS 
DUE  TO  THE  ADDITION  OF  DILUTE  POLYMER  SOLUTIONS 

J.R.  Bertshy  and  F.H.  Abernathy. 

Harvard  University,  U.S.A. 


Summary 

Two-dimensional  boundary  layer  flows  of  dilute  solutions  of  Polyox,  WSR-301,  have  been 
investigated  on  a free  surface  water  table.  The  natural  disturbance  level  was  low  enough  so 
imdisturbed  water  flows  were  laminar.  The  flows  had  wall  shear  velocities  up  to  4.2  cm/sec 
and  were  at  concentrations  up  to  32ppm  as  determined  by  testing  sanq>les  in  a vortex  viscometer. 
Observations  were  made  of  the  free  surface,  injected  dye,  streak  spacing  and  the  time-dependent 
velocity  in  the  mean  flow  direction  as  a function  of  polymer  concentration,  wall  strain  rate  and 
disturbance.  Nearly  continuous  and  very  localized  velocities  were  measured  using  a laser- 
Doppler  anemometer  and  a novel  method  of  selective  seeding.  Velocity  spectra  were  also 
obtained. 

Undisturbed  polymer  flows  at  sufficient  concentration  and  strain  rate  exhibited  velocity 
fluctuations  distinct  from  turbulence.  We  have  termed  these  as  polymer  Induced  fluctuations 
and  believe  that  they  are  most  likely  responsible  for  the  phenomenon  of  so  called  "early 
turbulence".  Studies  of  turbulent  polymer  flows  (disturbed  artificially)  indicated  significant 
drag  reduction  and  were  all  consistent  with  the  view  that  hairpin  eddies  are  the  major  mixing 
mechanism  in  a turbulent  boundary  layer,  and  that  these  eddies  are  sui^ressed  by  the  addition 
of  polymers.  When  strain  rate  and/or  concentration  were  increased  somewhat,  polymer 
fluctuations  were  found  to  coexist  with  the  turbulent  motions.  At  still  higher  strain  rates  and/or 
concentrations,  the  polymer  fluctuations  dominated  and  the  turbulent  flow  reached  suppressed 
condition.  No  distinction  between  artificially  disturbed  and  undisturbed  flows  at  the  same 
conditions  could  be  observed. 


Held  at  St.  John'i  Collegt,  Cambridge,  England. 

Organiied  and  iponaored  by  BHRA  Fluid  Engineering,  Cranfiatd,  Bedford,  MK43  OAJ. 
C BHRA  Fluid  Engineering. 
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NOMENCLATURE 


u ~ time  dependent  velocity  In  the  mean  flow  direction 
u'"*"  • rms  fluctuating  velocity  divided  by  u^ 

u"*”  “ mean  velocity  divided  by  u^ 

u^  “ wall  shear  velocity 

y"*"  • y u^/v  , y = distance  from  the  wall  v=solvent  viscosity  (water) 

“ Xu^/v  > “ mean  streak  spacing 


INTRODUCTTQN 

The  turbulent  drag  reducing  capability  of  dilute  polymer  solutions  has  been  known 
for  a number  of  years.  The  two  review  articles  by  Lumley  (ref.  1,2)  and  the  extensive 
surveys  of  Hoyt  (ref.  3)  and  Virk  (ref.  4)  review  and  place  in  perspective  the  known 
aspects  of  dilute  polymer  flows.  Yet,  mechanisms  by  which  the  turbulent  mixing  is 
reduced  are  still  poorly  understood.  For  example,  Rudd  (ref.  5),  Logan  (ref.  6)  and 
Kumor  and  Sylvester  (ref.  7)  have  each  reported  that  the^magnitude  of  the  fluctuation 
component  of  the  velocity  in  the  mean  flow  direction,  u'  , in  a turbulent  boundary 
layer  very  near  the  wall  increases  with  the  addition  of  polymers.  This  has  been 
interpreted  as  the  polymers  working  primarily  to  thicken  the  viscous  sublayer.  On  the 
other  hand,  Reischman  and  Tlederman  (ref.  8)  have  reported  a decrease  in  the  maximum 
intensity  of  u'"*"  near  the  wall  and  (along  with  ref.  9)  fo^r.d  a broad  peak  from  y‘^=  25 
and  beyond.  They  have  interpreted  this  as  an  indication  that  polymers  in  drag  reducing 
flows  primarily  alter  the  buffer  zone  (ref.  A).^.^5iey  reported  finding  no  indication 
of  a thickening  of  the  sublayer.  The  presepj^lnvestigation  was  designed  to  explore  the 
cause  of  these  observed  velocity  fluctuations  In  dilute  polymer  boundary  layers  and  to 
determine  in  what  ways  they  differed  from  those  observed  with  Newtonian  fluids. 

At  the  beginning  of  the  present  investigation,  we  held  the  view  that  the  major 
turbulent  mixing  mechanism  near  a solid  boundary  was  due  to  three-dimensional  flows 
induced  by  hairpin  eddies  aligned  in  the  flow  direction.  The  addition  of  dissolved 
pol'^nners  to  the  flow  should  inhibit  the  formation  and  hence  the  mixing  arising  from 
such  vortex  pairs.  At  high  enough  polymer  concentrations,  the  flow  near  the  wall 
should  be  essentially  laminar.  The  results  of  our  experiments  required  only  minor 
modification  of  this  initial  view  of  the  turbulent  boundary  layer. 

The  experiments  were  designed  to  investigate  boundary  layer  flows,  at  high  enough 
Reynolds  numbers  to  achieve  turbulence  if  artificially  disturbed,  but  at  low  enough 
free  stream  disturbance  levels  to  enable  a laminar  boundary  layer  flow  to  exist  in  the 
absence  of  such  disturbances.  This  enabled  us  to  observe  what,  if  any,  velocity 
fluctuations  might  arise  from  the  polymers  themselves. 


A free  surface  water  table  was  selected  as  the  optimal  flow  apparatus.  In 
addition  to  fulfilling  the  prerequisites  mentioned  above,  this  apparatus  is  well  suited 
for  optical  probing,  since  all  surfaces  are  planar  and  the  flow  is  easily  accessible. 
The  water  table  also  has  the  tremendously  useful  benefit  of  providing  information 
concerning  where  and  when  transition  to  turbulence  or  some  other  unsteady  flow  has 
occurred  by  visual  observation  of  the  free  surface. 


EXPERIMENTAL  FLOW  ARRANGEMENT  AND  EQUIPMENT 


A diagram  of  the  flow  arrangement  is  shown  in  Fig.  1.  Basically,  filtered  tap 
water  is  pumped  from  a sump  reservoir  through  an  overflow  standpipe  to  a water  table 
inclined  to  the  horizontal.  It  then  returns  to  the  sump  via  a weighing  tank.  Polymer 
solutions  were  injected  for  short  periods.  It  has  been  assumed  that  the  polymers  have 
been  severely  degraded  before  they  return  for  a second  time  to  the  water  table, 
especially  since  a large  depth  filter  (5p  nominal  pore  size)  was  interposed  in  the 
flow  circuit.  In  addition,  the  sump  reservoir  is  so  large  (approximately  17,000{.) 
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that  effects  of  polymers  returning  for  a second  pass  were  neither  expected  nor 
observed.  Concentrated  (200-500  ppm)  polymer  solutions  were  injected  by  gravity  feed 
from  constant  pressure  tanks  (two  2082.  drums)  into  the  piping  over  38  diameters  from 
the  water  table.  Thus,  the  final  dilution  to  the  test  concentration  was  accomplished 
by  turbulent  mixing  in  this  pipe  (10  cm  dla.). 

The  water  table  has  a deeper  inlet  section  equipped  with  large  cell  honeycomb  flow 
stralghteners  and  surface  wave  calming  devices.  There  is  a smooth  transition  to  the 
glass  working  section  which  is  1.22m  x 1.83m.  Typically,  the  flow  was  run  1.9  to 
3.5  mm  deep  and  the  table's  Inclination  was  from  1.6°  to  A. 2°.  All  undisturbed  water 
flows  presented  in  this  paper  were  laminar.  To  trip  this  flow  to  turbulence,  a bead- 
chain  of  2.63  mm  beads  was  inserted  in  the  flow  near  the  leading  edge  of  the  glass 
working  section. 

Velocity  measurements  were  made  using  a LDA  manufactured  by  Thermo  Systems,  Inc. 
(TSI)  oriented  as  shown  in  Fig.  2.  The  laser  used  had  a 1 mW  polarized  output. 

Signal  processing  was  accomplished  by  a TSI  model  1090  tracker.  It  is  normally 
difficult  to  resolve  the  local  velocity  in  a thin  shear  flow  using  LDA  techniques  since 
the  scattering  volume  of  the  intersecting  laser  beams  is  relatively  large  (here,  0.15mm 
dia.  by  1.7mm  length).  We  have  overcome  this  difficulty  by  introducing  light  scatter- 
ing contaminants  to  the  flow  only  at  a selected  distance  from  the  wall.  This  is 
accomplished  by  generating  bubbles  by  electrolysis  from  a fine  wire  held  parallel  to 
the  wall  and  normal  to  the  flow.  Most  other  scattering  particles  are  filtered  out  of 
the  flows.  Hence,  using  a forward  scatter,  fringe  mode  LDA  with  the  long  axis  of  the 
scattering  volume  normal  to  the  wall,  we  can  obtain  laminar  and  turbulent  velocity 
profiles  relatively  rapidly  by  moving  the  bubble  generating  wire  normal  to  the  wall 
(cf.  Fig.  2). 

The  bubble  wire  used  was  80%Pt20%Ir , 0.0025  cm  in  diameter  and  5 cm  long.  It  was 
supported  across  two  insulated,  L-shaped,  stainless  steel  tubes  which  entered  the  flow 
downstream  of  the  scattering  volume.  The  wire  was  always  positioned  at  least  300 
diameters  upstream  of  the  scattering  volume  so  that  its  wake  would  not  be  detected 
(ref.  10).  The  flow  was  always  such  that  the  local  Reynolds  number  of  the  wire  was 
less  than  25,  so  instability  of  the  wire  wake  was  never  a problem.  By  adding  a salt 
to  the  water  (Na2S0^,  anhydrous)  we  were  able  to  maintain  a sufficiently  high  bubble 
production  rate  so  that  the  measured  velocity  was  practically  continuous  (from  1000  to 
6000  velocity  samples/sec.)  while  the  operating  voltage  was  reasonably  safe  (<300V).* 

To  obtain  mean  velocities  and  rms  turbulent  intensities,  a signal  analysis  box 
developed  at  the  lab  was  used.  It  measured  the  mean  voltage  of  the  tracker's  output 
by  a continuously  integrating  op-amp  with  an  accuracy  of  0.1%.  This  voltage  was 
electronically  subtracted  from  the  tracker  output  and  then  fed  to  a Burr-Brown  4128  rms 
module  which  determined  true  rms  of  this  difference.  Its  accuracy  is  about  0.3%  from 
0.1  Hz  to  20  kHz.  Readings  were  taken  using  a 3^  digit  Data  Precision  digital  multi- 
meter. Due  to  the  one-pass  nature  of  the  experiments,  short  time  constants  were  used 
ana  the  measurements  were  considerably  more  uncertain  than  the  potential  accuracy  of 
this  instrumentation. 

Velocity  spectra  were  taken  using  a Tektronix  7L5  Spectrum  Analyzer.  It  operates 
digitally  and  was  set  to  filter  with  a 10  Hz  bandwidth.  Since  it  could  not  accept  a 
dc  signal  from  a high  Impedance  source,  a capacitor  coupled  op  amp  of  gain  1 was  inter- 
posed between  the  tracker  and  the  spectrum  analyzer. 

Depth  measurements  were  taken  using  a needle,  point  gauge  accurate  to  0.001  cm  for 
flows  with  a smooth  free  surface.  The  water  table's  inclination  was  monitored  by  a 
254  cm  long  liquid  level  located  along  the  table  side.  A thermometer  in  the  stilling 
section  read  temperature  with  an  accuracy  of  0.1°C. 


* Cambridge  City  water  contains  typically  100-150  ppm  of  dissolved  salts.  The 
addition  of  100-150  ppm  of  Na2S0^,  to  increase  the  electriaal  conductivity,  did  not 
effect  the  occurrence  of  polymer  induced  fluctuations  as  determined  by  viewing  the  free 
surface.  While  no  quantitative  measurements  of  the  dependence  of  polymer  effects  on 
salt  concentration  were  made,  we  have  no  evidence  to  suggest  the  results  were 
influenced  by  the  salts. 
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Polymer  Preparation  and  the  Rheological  Test  of  Concentration 

Batches  of  O.IZ  solution  (1511)  of  PEO,  WSR-301,  were  prepared  by  sprinkling  the 
dry  powder  on  the  surface  and  hand  mixing  gently  with  a large  paddle.  This  solution 
was  allowed  to  sit  for  two  days  before  using.  The  dry  polymer  was  not  dispersed  In  an 
alcohol  slurry  since  the  flow  over  rhe  water  table  Is  shallow  enough  that  It  might  be 
sensitive  to  surface  tension  changes.  Surface  tension  measurements  Indicated  that 
polymer  solutions  of  much  higher  concentrations  than  those  used  In  the  water  table 
exhibited  negligible  deviation  from  that  of  water. 

The  polymer  concentration  of  a given  flow  was  determined  by  withdrawing  a sample 
during  each  run  and  observing  its  behavior  in  a large  scale  vortical  flow  similar  to 
that  of  Balakrlshnan  and  Gordon  (ref,  11).  They  have  demonstrated  that  measuring  the 
minimum  polymer  concentration  needed  to  achieve  "vortex  Inhibition"  Is  a more  sensitive 
test  for  degradation  than  intrinsic  viscosity  measurements.  So,  we  used  a modified 
form  of  their  bathtub  vortex  experiment  as  a rheological  test  for  polymer  concentration. 
Briefly  th^  procedure  was  as  follows:  201  of  the  fluid  to  be  tested  was  placed  In  a 
cylindrical  container  (29.2  cm  inside  dia.)  with  a central  drain.  The  container  was 
rotated  until  solid-body  rotation  was  achieved  (50  rpm)  and  then  quickly  stopped. 
Draining  was  commenced  immediately  and  the  vortex  was  allowed  to  form  under  the 
maximum  draining  rate.  The  drain  was  connected  to  a long  flexible  tube  whose  outlet 
elevation  could  be  freely  varied. 

Since  nearly  all  the  flows  we  ran  over  the  water  table  had  concentrations  below 
"vortex  inhibition"  (25  ppm  for  WSR-301)  the  air-liquid  Interface  reached  to  the  drain 
and  there  was  two-phase  flow  in  the  exhaust  tube.  As  soon  as  the  vortex  was  securely 
established,  the  end  of  the  exhaust  tube  was  raised  until  two-phase  flow  no  longer 
occurred  in  it.  This  elevation  was  then  compared  to  a calibration  curve  made  using 
dilute  solutions  made  from  a single  batch  of  0.1%  stock  solution.  Thus,  we  were  able 
to  judge  the  concentration  compared  with  this  single  stock  solution  to  about  2 ppm 
independently  of  how  poorly  prepared  or  mishandled  the  sample  had  been.  Typically, 
samples  after  passing  over  the  water  table  behaved  in  the  vortex  as  if  their 
concentration  had  been  25%  to  50%  less  than  that  injected. 

The  effective  concentration  of  polymer  determined  by  the  vortex  viscometer  was  in 
agreement  with  the  water  table  flow  observations.  All  data  is  presented  with  the 
concentration  as  determined  in  this  fashion.  In  the  near  future,  we  expect  to  under- 
take a detailed  study  of  this  vortex  flow. 

Velocity  Data 

In  so  far  as  time  allowed,  complete  profiles  of  u-velocity  and  turbulence 
intensity,  were  measured  for  each  polymer  run  (45  minutes  or  less).  Since  the  LDA 
would  occasionally  lose  track  of  the  Doppler  signal,  only  short  time  averaging 
(about  8 sec.)  was  possible  to  ensure  quick  recovery.  In  addition,  bubble  production 
was  reduced  in  polymer  flows  so  constant  attention  was  necessary  to  obtain  data. 

All  velocity  and  depth  measurements  were  taken  at  a station  146  cm  downstream  of 
the  transition  from  the  stilling  section  to  the  test  section.  Since  the  flows  were 
typically  0.2  cm  to  0.3  cm  deep,  this  corresponds  to  a test  location  between  730  and 
500  flow  depths  downstream  of  the  transition. 

Wall  shear  stress  and  hence  u^  was  inferred  from  the  depth  measurement  and  the 
table  inclination  and  assuming  flow  equilibrium.  The  flow  accelerates  as  it  progresses 
down  the  table  and  wall  shear  stress  accordingly  increases.  A turbulent  water  flow 
comes  to  an  equilibrium  depth  (and  hence  velocity)  within  60  cm  cn  the  flat  test 
section,  while  laminar  water  flow  has  not  quite  come  to  equilibrium  even  at  the 
measuring  station  146  cm  downstream.  Estimates  of  the  acceleration  of  laminar  water 
flows  similar  to  those  in  which  polymers  were  tested  indicate  that  at  most,  might 
be  decreased  by  4%  from  the  value  used.  Flows  indicating  any  fluctuation  and  hence 
an  increase  in  momentum  exchange  from  laminar  water  flows  were  certainly  much  closer 
to  equilibrium.  Uncertainty  in  depth  and  angle  measurements  lead  to  an  uncertainty 

in  u of  + 0.02  cm/sec  for  laminar  flows  and  -I-  0.06  cm/sec  for  turbulent  flows. 
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The  position  of  the  bubble  wire  relative  to  the  wall  was  also  determined  by  the 
depth  and  angle  measurements.  Average  velocity  data  from  a laminar  flow  was  plotted 
with  relative  uncertainty  In  position  of  only  + 0.001  cm.  The  wall  shear  stress 
measure  (assuming  solvent  kinematic  viscosity)  was  then  used  as  an  Indication  of  the 
slope  this  velocity  data  had  to  have  at  the  wall.  The  two  were  then  compared  and  an 
absolute  measure  of  distance  from  the  wall  to  the  bubble  wire  was  determined  with  an 
uncertainty  of  approximately  + 0.005  cm. 

All  data  Is  presented  using  the  solvent  kinematic  viscosity  as  determined  by  the 
temperature.  There  Is  evidence  that  this  may  not  be  the  best  method  to  correlate  the 
data;  however,  most  previous  investigators  have  presented  their  data  using  solvent 
viscosity.  We  do  so  to  facilitate  comparison,  but  point  out  that  Relschman  and 
Tlederman  (ref.  8)  did  not  follow  this  practice  and  allowances  should  be  made 
accordingly. 

One  further  flow  visualization  experiment  was  performed  to  obtain  the  streak 
spacing.  Grains  of  Na-SO.  were  sprinkled  Into  the  flow.  Since  the  salt  Is  denser 
than  water,  it  tended  ^to^travel  downstream  along  the  bottom  glass  wall  of  the  water 
table.  The  grains  are  small  enough,  however,  to  be  lifted  by  turbulent  motions  and 
streaks  were  readily  observed.  Since  the  terminal  velocity  of  these  salt  grains  In 
still  water  Is  about  1 or  3 cm/sec.  It  Is  likely  that  the  turbulent  motions  must 
occasionally  exceed  this  velocity  for  streaks  to  be  detected  In  this  manner. 

EXPERIMENTAL  RESULTS 


The  philosophy  used  throughout  this  Investigation  has  been  to  compare  the  dilute 
polymer  flows  with  those  of  pure  water  at  comparable  strain  rates  and  depths  of  flow. 
Following  such  a program  required  a careful  investigation  of  the  flow  of  water, 
measuring  all  of  the  variables  with  the  same  techniques  to  be  used  with  dilute  polymer 
flow. 


The  dimensionless  mean  velocity  U is  plotted  versus  dimensionless  distance  away 
from  the  wall  y*^  for  both  laminar  and  turbulent  flow  of  water  in  Fig.  3.  The  measured 
results  for  the  turbulent  boundary  layer  flow  of  water  show  very  close  agreement  with 
the  "law  of  the  wall"  results  for  Newtonian  fluids  from  both  boundary  layer  and  pipe 
flow.  In  fact,  the  water  table  turbulence  data  Is  nicely  bracketed  by  the  extremes 
commonly  quoted  in  the  literature  for  this  universal  profile  (ref.  12).  Departure  from 
the  universal  velocity  profile  naturally  does  occur  near  the  free  surface  since  the 


slope  ^ must  go  to  zero  there.  The  laminar  profile  is  a reasonable  approximation  to 

the  calculated  asymptotic  l^lt.  The  nondimens ional  rms  value  of  the  fluctuating 
component  of  u velocity,  u'  Is  plotted  In  Fig.  4 for  turbulent  water  flow_^along  with 
the  laminar  value,  which  results  from  the  noise  level  of  the  LDA.  (All  u'^  values 
presented  have  been  corrected  for  this  level.)  The  general  shape  and  the  location  of 
the  peak  of  u'"*"  at  y"*"  “^15  conforms  with  the  published  measurements  (ref.  12);  however, 
the  overall  scale  of  u'  is  lower  for  the  water  table  flow.  This  should  be  expected 
because  the  water  table  flow  extends  only  to  y*^  on  the  order  of  100,  while  pipe  flow 
and  boundary  layer  flow  generally  have  a log  layer  flow  several  orders  of  magnitude 
larger.  The  limited  log  lay^r  of  water  table  flow  does  not  effect  the  mean  profile 
(Fig.  3)  nor  Che  shape  of  u'  , It  only  decreases  the  latter's  scale. 


The  mean  and  fluctuating  velocity  profiles  for  water  (Figs.  3 and  4)  were 
obtained  using  the  selective  seeding  LDA  technique.  Nearly  continuous  velocity 
measurements  can  be  frequency  analyzed,  and  Fig.  5 shows  the  oscilloscope  output  from 
the  spectrum  analyzer  at  y***  - 14,  48  and  58  for  a water  turbulent  boundary  layer. 

These  three  values  of  y are  In  the  neighborhood  of  those  that  will  be  presented  for 
polymer  flows.  The  fluctuations  in  spectral  Intensity  are  due  to  the  manner  In  which 
the  analyzer  operates.  Its  10  Hz  pass  band  slowly  and  automatically  sweeps  across  the 
frequency  range  of  0 to  500  Hz.  When  centered  at  a particular  frequency,  the  signal 
then  entering  Che  Instrument  Is  digitally  analyzed.  Successive  averaging  of  such 
spectra  would  naturally  decrease  Che  fluctuations  In  Fig.  5,  resulting  In  a smooth 
curve  over  the  0 to  500  Hz  range.  The  three  spectral  traces  In  Fig.  5 show  no  major 
variations  from  y ■ 14  to  y ■ 54,  suggesting  Chat  the  same  disturbances  are  active 
throughout  the  buffer  layer.  A large  fraction  of  the  fluctuation  is  at  low  frequency. 


below  75  Hz.  Assuming  these  fluctuations  are  convected  with  the  local  velocity 
suggests  that  their  wavelength  Is  greater  than  the  thickness  of  the  layer,  which  Is 
verified  by  visual  observations  of  the  free  surface. 

The  object  of  the  experiments  was  to  determine  how  polymer  concentrations  and  wall 
strain  rate  modify  laminar  flows  and  then  to  induce  transition  to  turbulence  and  again 
measure  Che  same  parameter  dependence.  Since  the  wall  shear  stress  as  of  turbulent  and 
laminar  flow  at  constant  mass  flow  and  table  inclination  differ  In  general,  a patch  of 
turbulent  flow  will  be  deeper  than  a laminar  region  and  therefore  easily  seen  by 
looking  at  the  free  surface  (Fig.  6a).  It  was  this  property  which  led  to  the  discovery 
by  Emmons  (ref.  13)  of  spots  on  a water  table.  Such  spots  can  be  reproduclbly  generated 
by  Injecting  a pulse  of  fluid  from  a hypodermic  Cube  upstream  of  the  entrance  contrac- 
tion. The  Injected  fluid  creates  a three-dimensional  disturbance  (mostly  a vortex  ring) 
which  triggers  a localized  transition  to  turbulence.  The  spanwlse  extent  of  the 
turbulence  Is  many  times  the  depth  of  Che  flow.  In  water  the  width  of  the  spot  along 
the  span  Increases  as  the  spot  Is  convected  downstream.  Recently  (ref.  14)  spots  In  an 
air  boundary  layer  have  been  studied  and  Che  velocity  profile  In  the  spot  away  from  the 
leading  and  trailing  edges  has  been  found  to  be  fully  turbulent. 

With  the  addition  of  polymers  to  the  flow,  spot  formation  Is  altered  dramatically 
(Fig.  6b).  In  contrast  to  a spot  In  pure  water  the  lateral  growth  rate  Is  greatly 
reduced  and  the  depth  of  a spot  Is  not  much  different  from  the  laminar  flow  surrounding 
It.  As  the  shear  stress  or  polymer  concentration  Increases,  fluctuations  on  the  free 
surface  of  the  otherwise  laminar  flow  appear  and  can  be  seen  in  the  background  flow  of 
Fig.  6c.  These  fluctuations  are  distinct  from  those  due  to  turbulence,  since  transi- 
tion to  turbulence  can  still  be  Induced  (Fig.  6c).  An  Increase  In  wall  shear  (above 
the  laminar  value)  In  pipe  flows  at  Reynolds  numbers  below  the  expected  turbulent 
transition  has  been  observed  (ref.  15)  and  this  was  called  early  turbulence.  It 
appears  likely  that  this  Increase  In  drag  can  be  attributed  to  the  velocity  fluctuations 
visible  in  Fig.  6c.  We  have  called  these  velocity  fluctuations  polymer  Induced 
fluctuations  because  they  differ  from  turbulence.  The  distinction  between  polymer 
Induced  velocity  fluctuations  In  laminar  flow  and  turbulence  fluctuations  Is  sharp  at 
moderate  concentrations  and  strain  rates  and  becomes  less  distinct  at  higher  concen- 
trations and  strain  rates.  Hence,  It  will  be  Important  to  distinguish  between  laminar- 
like  properties  and  turbulent-llke  properties.  For  laminar  flow: 

(1)  The  free  surface  Is  smooth  and  free  of  waves  (background  of  Fig.  6a,  away  from  the 
spots) . 

(2)  The  u-veloclty  trace  Is  relatively  smooth;  Its  spectrum  is  broadband,  and  u’'*'  Is 
low  and  essentially  Independent  of  y*^  (Fig.  4). 

(3)  Dye  Introduced  at  the  wall  Is  convected  along  the  wall  with  little  lateral 

spread;  salts  Introduced  in  the  flow  show  no  evidence  of  streaks. 

(4)  By  Introducing  disturbances,  either  by  a trip  cord  or  by  Injecting  fluid,  an 

Instability  can  be  generated  which  causes  transition  to  turbulence. 

(5)  The  mean  velocity  profile  can  be  distinguished  from  the  turbulent  profile 
(Fig.  3 for  water). 

For  turbulent  flow; 

(1)  The  free  surface  Is  wavy  with  many  wavelengths  greater  than  depth  (the  spot  In 
Fig.  6a). 

(2)  The  u-velocity  has  a large  fluctuating  component  (Fig.  4)  and  spectra  (Fig.  5) 
have  large  low  frequency  components  and  are  essentially  Independent  of  y"^. 

(3)  Dye  introduced  at  the  wall  la  rapidly  dispersed;  s^lts  introduced  at  the  wall 
give  evidence  of  streaks  with  a spanwlse  spacing  (X  ) of  116  ± 15. 

(4)  Introducing  disturbances  has  no  noticeable  effects  on  the  flow  downstream  of  the 
disturbance. 


01-6 


(5)  The  mean  velocity  profile  follows  the  universal  velocity  distribution  for 

turbulent  flow  (Fig.  3). 

Laminar  Flow  with  Polymers 

Laminar  flows  with  dissolved  polymers  exhibit  fluctuations  on  the  free  surface 
only  for  certain  values  of  concentration  and  wall  strain  rate  (dU/dy),  as  shown  In 
Fig.  7.  No  attempt  was  made  to  define  the  region  precisely  as  It  was  enough  for  this 
Initial  Investigation  to  determine  the  general  boundaries  of  the  region.  The  form 
that  this  boundary  curve  assumes  Is  consistent  with  observations  for  the  onset  of 
drag  reduction.  Berman  and  George  (ref.  16)  suggest  a time  scale  criterion  which  for 
WSR-301  predicts  an  onset  strain  rate  of  about  320  sec”^*,  certainly  less  than  the 
strain  at  the  knee  720  sec"^  of  Fig.  7.  While  drag  reduction  has  been  observed  near 
the  limit  of  Infinite  dilution,  the  onset  strain  rate  Increases  as  concentration 
decreases  (ref.  17).  This  too  Is  In  accord  with  the  trends  of  Fig.  7,  suggesting  that 
polymer  fluctuations  are  closely  related  to  the  mechanisms  of  polymer  drag  reduction. 

Graphs  of  u'"*"  vs  y"*"  for  polymer  fluctuations  at  several  strain  rates  and  concentra- 
tions are  shown  In  Fig.  8.  At  first  glance,  they  are  similar  to  those  for  water 
turbulence  In  Fig.  4,  For  all  cases  except  for  the  highest  concentration  and  strain 
rate,  u*'*’  peaks  at  y =15  and  then  falls  off  to  a much  lower  value  at  higher  y . The 
strong  difference  between  the  polymer  fluctuation  and  water  turbulence  can  be  seen  In 
the  spectra  shown  In  Fig.  9 which  correspond  to  conditions  2 In  Fig.  8.  These  spectra 
show  the  features  common  to  all  cases  of  polymer  fluctuations,  namely  strong  spectral 
variation  with  y"^  while  the  low  and  high  frequencies  are  suppressed  compared  to  water 
turbulence.  Compared  to  laminar  flow  of  water,  there  Is  the  addition  of  fluctuating 
components  in  the  50  Hz  to  150  Hz  range.  With  these  fluctuations  confined  in  a region 
around  y » 15,  it  is  natural  to  assume  they  are  convected  with  the  local  mean  velocity 
implying  a wavelength  on  the  order  of  the  layer  thickness  or  larger.  An  Inspection  of 
the  free  surface  (cf.  Fig.  6c)  shows  wavelengths  on  the  order  of  the  layer  depth.  It 
Is  a characteristic  of  laminar  layers  with  polymer  fluctuations  to  have  the  more  Intense 
fluctuations  confined  to  a region  near  the  wall.  As  either  the  concentration  or  strain 
rate  Increases,  the  vertical  extent  increases,  out  to  y"*"  =35  at  the  highest  strain  rate 
and  concentrations  tested.  In  all  cases  the  polymer  fluctuations  as  detected  In  u''*' 
were  In  the  frequency  range  of  50  to  150  Hz  and^could  be  defected  only  near  the  wall 
(cf.  Fig.  10).  In  addition,  it  appears  that  u'  at  large  y'*’  tends  to  Increase  with 
concentration. 

The  mean  velocity  profiles  u"*"  versus  y"*"  were  quite  sensitive  to  polymer  concentra- 
tion. The  velocity  curves  1 and  2 in  Fig.  10  correspond  to  conditions  1 and  3 on 
Fig.  8.  At  lower  concentration,  for  example  13  ppm,  where  polymer  fluctuations  can 
first  be  detected,  the  measured  profiles  fall  between  water  laminar  and  curves  1 and  2. 
If  the  mean  velocity  curves  had  been  corrected  for  the  kinematic  viscosity's  dependence 
on  polymer  concentration  and  strain  rate,  then  the  profiles  would  have  fallen  closer 
to  the  water  laminar  curve.  The  measured  mean  velocity  profiles  for  polymer  fluctxia- 
tlons  are  similar  to  those  of  the  laminar  flow  of  water,  with  an  apparent  higher 
viscosity.  This  holds  true  up  to  effective  polymer  concentrations  of  25  to  32  ppm 
where  It  Is  Impossible  to  distinguish  between  polymer  fluctuations  and  polymer  turbul- 
ence If  the  strain  rate  Is  high  enough  for  polymer  effects  to  be  observed  (cf.  Fig.  7). 

For  flow  described  as  polymer  fluctuations,  attempts  were  made  to  find  evidence 
of  streaks  In  the  flow  direction  but  none  were  observed.  For  water  turbulence,  streaks 
were  observed,  and  their  mean  spacing  In  dimensionless  wall  variables  was  116  + 15,  In 
agreement  with  measurements  In  air  and  water  of  many  others  using  many  different 
techniques. 

Turbulent  Flows  with  Polymers 

The  same  type  of  experiments  and  measurements  already  reported  for  the  polymer 
fluctuations  were  repeated  for  turbulent  flows.  It  was  necessary  to  Introduce  a 


*Berman  and  George  (ref.  16)  used  WSR-N80  and  determined  onset  for  glycerine-water 
polymer  solutions.  We  have  adjusted  their  value  to  account  for  the  higher  molecular 
weight  and  Intrinsic  viscosity  of  WSR-301  to  arrive  at  this  figure  for  onset. 
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larger  disturbance  to  cause  transition  to  turbulence  at  higher  polymer  concentrations 
than  it  was  for  pure  water  transition.  The  larger  disturbances  were  either  longer  time 
pulses  of  fluid  to  cause  transition  in  spots  (cf.  Fig.  6b)  or  larger  diameter  rods  that 
spanned  the  flow  at  the  entrance  station.  Hansen  et  al  (ref.  18)  on  the  basis  of  a 
linear  stability  analysis  of  pipe  flow,  using  a single  time  scale  viscosity  model  for 
dilute  polymers,  calculated  that  the  growth  rate  of  disturbances  would  be  altered  by 
the  polymers.  From  our  observations  of  the  water  table  flows,  we  can  report  that 
larger  amplitude  disturbances  are  definitely  necessary  for  transition,  and  that  the 
spanwlse  growth  of  turbulent  spots  is  also  greatly  reduced.  In  fact,  three-dimensional 
spanwlse  variations  in  disturbances  appear  to  be  required  to  insure  reproducible 
turbulent  flow.  A bead-chain  (in  fact  an  aluminium  electric  light  pull  chain)  was  used 
which  completely  spanned  the  flow  at  the  entrance  station.  The  bead-chain  seemed 
satisfactory  for  all  flows.  Observations  were  made  of  the  free  surface  with  and  with- 
out the  chain.  Additional  disturbances  such  as  cylinder  2 or  3 cm  in  dla.  placed  in 
the  flow  downstream  of  the  chain  resulted  in  no  observable  effects  for  more  than  a few 
diameters.  It  was  therefore  concluded  that  turbulent  conditions  with  polymers  had  been 
achieved,  and  it  is  this  condition  which  was  used  for  the  reported  flow  measurements. 

The  visual  appearance  of  free  surface  is  different  for  polymer  turbulence  and 
polymer  fluctuations  at  concentrations  below  25  ppm  as  can  be  seen  in  Figs.  6b  and  6c. 
The  amplitude  and  wavelength  of  the  free  surface  fluctuations  ace  much  larger  for 
turbulence  than  for  polymer  fluctuations,  as  the  spectrum  of  u*'*'  reveals.  In  addition, 
the  turbulent  patch  is  deeper  than  the  surrounding  flow  since  the  wall  shear  stress 
under  the  patch  is  greater. 

The  measurements  of  u'^  vs  y^  for  several  polymer  concentrations  and  strain  rates 
are  plotted  in  Fig.  11  along  with  that  of  pure  water  for  reference.  At  the  lowest 
concentration  there  is  an  increase  in  the  magnitude  of  u''*’  over  that  of  pure  water, 
even  though  the  spectral  distribution  (Fig.  12)  is  quite  similar  to  that  of  pure  water. 
At  this  concentration,  no  polymer  fluctuations  were  detected  in  the  laminar  flow. 

Curve  2 in  Fig.  11  shows  a double  peak  in  intensity,  and  its  u-veloclty  spectra  at  four 
different  values  of  y'*’  are  shown  in  Fig.  13.  The  first  peak  near  y"*”  « 15  appears  to  be 
associated  with  polymer  fluctuations  both  on  the  basis  of  its  y'*'  location  and  the  u- 
veloclty  spectrum.  The  increase  in  spectral  intensity  at  lower  frequency  is  definitely 
present.  The  second  peak  near  y"^  of  30,  on  the  basis  of  the  spectra  at  y"*"  “ 32  and  41, 
is  probably  associated  with  the  turbulence  process. 

From  the  wide  range  of  the  u'"*"  versus  y"*"  data  as  a function  of  polymer  concentra- 
tion and  strain  rate,  it  is  apparent  that  a simple  characterization  of  the  fluctuation 
dependence  is  not  possible.  Our  results  for  a two-dimension  mean  flow  indicate  a 
strong  dependence  on  wall  strain  rate  over  the  entire  range  of  polymer  concejjitratlon 
from  low  to  intermediate,  and  to  high.  Curve  1 of  Fig.  11  indicates  that  u'  is 
augmented  over  the  solvent  values  at  low  concentrations;  curves  2,  3 and  4 are  all  for 
a moderate  concentration  of~16  ppm  and  show  the  richness  of  the  strain  rate  dependence. 
The  lower  strain  rate  (curve  3;  u ■ 3.25  cm/sec)  has  u''*'  below  the  solvent  value  for 
all  y'*'.  The  u velocity  spectra  at  the  next  highest  strain  rate  (u  « 3.52  cm/sec) 
in  Fig.  13  shows  evidence  of  strong  polymer  fluctuations  near  y***  *^15,  the  first  peak 
in  the  u'  intensity.  Increasing  u to  3.85  cm/sec  lowers  the  intensity  of  u'"^  every- 
where, the  spectrum  near  the  peak  shown  in  Fig.  14  being  suggestive  of  polymer 
fluctuations.  Curves  5 and  6 of  Fig.  11  are  representative  of  what  we  have  called 
asymptotic  conditions  or  suppressed  turbulence.  The  effective  concentrations  and 
strain  rates  are  so  large  that  no  distinction  can  be  made  between  the  flow  with  and 
without  the  bead-chain  disturbance.  Spectra  for  curve  5 of  Fig.  11  are  shown  in 
Fig.  15.  The  flow  is  inactive  except  for  the  fluctuations  near  the  wall  which  which 
are  suggestive  of  polymer  fluctuations. 

The  turbulent  mean  velocity  profiles  for  concentrations  below  21  ppm  are  plotted 
in  Fig.  16.  The  profiles  for  the  asymptotic  cases  are  the  solid  circles  and  diamonds 
on  Fig.  10  which  are  plotted  with  other  polymer  fluctuation  flows.  The  asymptotic 
cases  fit  into  an  ordered  sequence,  as  concentrations  and  strain  rates  are  increased 
when  viewed  as  laminar  flows  with  polymer  fluctuations.  VHien  viewed  from  the  water 
turbulent  limit,  the  mean  velocity  profiles  increase  their  difference  from  the 
Newtonian  turbulence  as  concentration  increases,  but  then  fall  back  as  the  asymptotic 
conditions  are  approached  and  polymer  fluctuations  play  a dominant  role. 
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The  mean  velocity  profiles  plotted  in  the  standard  wall  coordinates  are  indicative 
of  the  effective  viscosity  in  the  different  ranges  in  the  flow.  In  the  sublayer, 
y < 8,  the  mean  velocity  profiles  for  laminar  and  turbulent  water  arc  the  same  and 
there  is  no  evidence  in  our  experiments  that  polymer  flows  with  either  polymer  induced 
fluctuations  or  turbulence  departs  significantly  from  this  solvent  behavior.  As  far  as 
the  sublayer  thickness  is  concerned,  our  results  are  in  agreement  with  those  of 
Relschman  and  Tiederman  (ref.  8);  however,  the  laminar  water  flow  on  the  table  departs 
from  ■ y*^  in  the  neighborhood  of  y***  ==8  (cf.  Fig.  3)  because  the  shear  stress  is  not 
constant  throughout  the  flow.  The  laminar  water  flow  mean  velocity  profile  represents 
what  is  expected  when  the  viscosity  is  constant  throughout  the  flow.  On  the  other  hand, 
the  turbulent  water  flow  profile  is  the  result  of  augmentation  of  molecular  viscosity 
by  the  turbulent  mixing  process.  Departures  from  this  profile  for  turbulent  polymer 
flows  imply  changes  in  the  mixing  process.  Following  this  line  of  reasoning.  Fig.  16 
shows  a departure  from  the  universal  profile  for  9 ppm  at  u «*  3.72  cm/sec  and  the 
maximum  departure  at  lA  ppm  and  >3.77  cm/sec.  At  this  latter  concentration, 
polymer  fluctuations  were  observed  at  u^  * 3.71  cm/sec, which  is  near  the  boundary  curve 
for  Fig.  7.  Increasing  polymer  concentration  and  strain  rate  appears  to  inhibit  normal 
turbulent  mixing;  however,  beyond  the  moderate  concentration  of  14  ppm  fluctuations 
in  the  u-veloclty  are  Induced  In  the  laminar  flow  resulting  in  an  effective  increased 
viscosity,  at  least  in  the  region  where  the  fluctuations  are  observed.  The  resulting 
turbulent  flows  at  concentrations  of  21  ppm  and  beyond  appear  not  only  to  consist  of 
suppressed  turbulent  mixing  but  also  have  augmented  viscosity  due  to  polymer  fluctua- 
tions. Those  turbulent  flows  with  the  lowest  p^k  values  of  u*'*'  appear  to  be  dominated 
by  polymer  Induced  fluctuations  observable  at  y near  the  peak  (Fig.  15)  rather  than  by 
normal  turbulent  mixing  processes  which  show  a similar  profile  of  u*  (cf.  Fig.  5 for 
pure  water  of  Fig.  12  for  9 ppm  of  polymer).  Increasing  the  polymer  concentration  from 
25  to  32  ppm  increases  the  maximum  of  u'  relative  to  its  value  at  15  ppm  and  shifts 
the  resulting  mean  turbulent  velocity  profile  back  toward  the  water  turbulent  curve. 

The  shift  is  due  to  an  increase  in  polymer  fluctuations  rather  than  due  to  Increased 
turbulent  mixing.  The  asymptotic  flow  situations  exhibit  more  laminar  features  than 
turbulent  ones.  S>pectra  for  this  flow  are  shown  in  Fig.  15.  The  fluctuations  are 
confined  to  the  wall  region,  outside  of  which  the  velocity  gives  the  appearance  of  a 
laminar  flow  with  possibly  presence  of  some  very  low  frequencies. 

Measurements  of  the  average  spanwlse  streak  spacing,  X'*',  were  made  of  the  polymer 
turbulent  flows  using  the  ^alt  crystal  technique.  A few  examples  illustrate  the 
results.  For  pure  water  X’’’  ■ 116  + 15;  at  concentration  of  9 ppm  and  u • 3.72  cm/sec, 
X"*"  " 161  + 115;  at  9 ppm  and  u » 4.21  cm/sec,  X"*"  • 200  + 18;  at  16  ppm  and  u^  • 3.52 
cm/sec  X+  ■ 190  + 20.  For  the  (lows  at  and  near  asymptotic  conditions  there  was  no 
evidence  of  streaks  using  this  technique. 


DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 


The  flows  investigated  were  two-dimensional  laminar  and  turbulent  shear  flows  on  a 
free  surface  water  table.  These  basic  flo|^s  were  characterized  in  detail.  Measurements 
of  U^,  u'^  and  u-velocity  spectra  versus  y have  been  reported  along  with  observations 
of  the  free  surface  and  measurements  of  spanwise  streak  spacing  for  turbulent  flows. 

The  laminar  flow  is  nearly  an  equilibrium  flow  and  in  good  agreement  with  the  calculated 
profile  assuming  the  body  force  is  equal  to  the  local  wall  shear  stress.  The  turbulent 
flow  mean  velocity  profile  is  in  excellent  agreement  with  the  universal  velocity  ^ 
profile  even  though  the  flow  extends  only  to  y'*’  on  the  order  of  100.  The  peak  of  u' 
occurs  at  y"*"  of  15  in  agreement  with  the  results  of  other  investigations  made  in  other 
turbulent  boundary  layers  and  in  pipes,  while  the  absence  of  the  outer  flows  appears 
to  just  reduce  the  overall  level  of  u'*^.  The  spanwlse  streak  spacing  is  in  agreement 
with  published  values  of  other  investigators.  The  absence  of  the  outer  flow  allowed 
the  near  wall  region  to  be  studied  without  the  extraneous  disturbances  Imposed  by  the 
outer  flow.  As  a result,  it  has  been  possible  to  establish  features  of  the  inner  flow 
which  have  not  been  revealed  in  other  investigations. 

At  modest  polymer  concentrations  and  strain  rates  large  enough  to  be  comparable  to 
the  time  scale  of  the  polymer,  effects  on  the  velocity  very  near  the  wall  become 
noticeable  in  an  otherwise  laminar  flow.  Th^  mean  intensity  of  polymer  Induced 
fluctuations  has  a spatial  distribution  in  y^  similar  to  what  one  finds  in  Newtonian 
turbulence,  though  the  intensity  is  lower  and  spectral  distribution  is  quite  different. 
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The  distinction  is  further  confirmed  by  velocity  time  traces  and  the  free  surface 
appearance.  No  evidence  of  streaks  in  the  flow  were  found  at  the  wall.  Transition  to 
a polymer-influenced  turbulence  would  be  triggered  by  a large  enough  three-dimensional 
disturbance  in  the  flow. 

The  polymer  Induced  velocity  fluctuations  are  clearly  strain  rate  dependent  as 
shown  in  Fig.  7.  Since  the  strain  rate  where  they  are  first  observed  is  of  the  same 
order  as  the  Inverse  of  the  calculated  response  time  for  the  polymer,  it  is  natural  to 
assume  these  fluctuations  are  associated  with  at  least  partial  expansion  of  the  polymer. 
Observations  of  an  Increase  in  laminar  wall  shear  stress  in  capillary  flows  at  similar 
maximum  strain  rates  have  been  ascribed  to  "early  turbulence"  (ref.  15).  We  have  also 
observed  an  Increase  in  wall  shear  stress  when  polymer  Induced  fluctuations  are  present; 
however,  this  flow  state  is  distinct  from  turbulence,  since  at  modest  polymer  concen- 
trations a separate  and  distinct  transition  can  be  Induced.  In  addition,  the  spatial 
variation  of  u-velocity  spectra  shows  that  polymer  Induced  fluctuations  are  primarily 
confined  to  the  region  near  the  wall  with  the  fluctuations  primarily  in  the  50  to  150 
Hz  range. 

Turbulent  flows  of  polymer  solutions  at  low  and  modest  concentrations  or  strain 
rates  did  exhibit  streaks  in  the  sublayer.  The  spacing  between  streaks  Increased  with 
concentration  and  strain  rate.  In  this  regime,  turbulent  fluctuations  were  similar 
throughout  the  layer,  and  the  spectra  were  similar  to  those  of  water.  By  comparing  the 
depth  of  the  polymer  flow  to  that  of  water  at  the  same  mass  flow  rate  and  table  incli- 
nation, drag  reduction  was  definitely  established.  In  addition,  the  mean  velocity 
profiles  showed  a departure  from  the  universal  velocity  profile,  suggesting  a change 
in  turbulent  mixing  or  effective  turbulent  viscosity.  Polymer  effects  on  the  mean 
turbulent  velocity  profile  and  on  wall  shear  stress  were  observed  at  polymer  concen- 
trations and  strain  rates  below  those  needed  to  cause  polymer  fluctuations.  Though  no 
attempt  was  made  to  define  an  onset  curve  for  polymer  effects  on  turbulence  similar  to 
the  one  for  polymer  fluctuations  in  laminar  flows,  the  experimental  data  suggest  that 
the  onset  line  lies  to  the  left  and  below  the  polymer  fluctuation  curve  of  Fig.  7. 

Such  a result  is  consistent  with  our  view  of  hairpin  eddies  being  the  major  mixing 
mechanism  in  a Newtonian  turbulent  boundary.  The  maximum  strain  rate  in  an  eddy 
aligned  with  the  flow  would  be  higher  than  that  at  the  wall,  since  the  vortex  tube 
forming  the  eddy  is  stretched  by  the  mean  flow.  In  addition,  the  maximum  strain  rate 
occurs  where  rotational  effects  are  small,  outside  the  vortex  core.  Hence,  expansion 
of  the  polymer  coll  could  occur  at  a lower  strain  rate  than  it  could  in  the  laminar 
shear  flow  at  the  wall  where  strain  and  rotation  are  in  balance  (ref.  2).  For  both  of 
these  reasons,  it  should  be  expected  that  polymers  would  begin  to  inhibit  turbulent 
mixing  at  lower  strain  rates  than  In  laminar  flows.  At  a fixed  concentration,  strain 
rate  dependence  should  also  be  expected.  Polymers  such  as  WSR-301  contain  molecules 
of  different  lengths,  and  the  largest  expand  at  lower  strain  rates  than  the  shorter 
molecules.  Independent  of  the  exact  structure  of  the  three-dimensional  mixing 
mechanism,  some  of  the  structures  will  Involve  higher  strain  rate  motions  than  others. 

It  is  these  more  intense  structures  which  will  be  Inhibited  at  the  lowest  strain  rates, 
leading  to  an  Increase  in  turbulent  streak  spacing  as  strain  rate  is  Increased.  Such 
reasoning  is  consistent  with  the  experimental  result  presented.  It  is  further  con- 
firmed by  comparison  of  the  u'  spectrum  across  the  flow  as  a function  of  concentration 
and  strain  rate.  Fig.  5 is  the  water  spectra  and  Fig.  12  is  for  9 ppm,  a concentration 
below  polymer  induced  fluctuations;  the  spectra  are  quite  similar,  though  drag  reduc- 
tion is  observed  and  the  mean  velocity  profile  has  been  Influenced  indicating  a change 
in  mixing.  Increasing  concentration  and  strain  rate.  Figs.  13  and  14  show  the  velocity 
fluctuations  substantially  reduced  except  for  y*^  < 20,  where  polymer  Induced  fluctua- 
tions are  found  in  laminar  flows.  With  further  increase  in  strain  rate  and  concentra- 
tion, an  asymptotic  regime  is  reached  where  most  of  the  turbulent  features  are 
suppressed  but  the  polymer  Induced  velocity  fluctuations  remain  (Fig.  15).  Under  these 
conditions,  streaks  were  never  observed.  These  suppressed  turbulent  flows  have  mostly 
laminar  characteristics,  and  it  was  difficult  to  distinguish  these  flows  from  those  of 
polymer  induced  fluctuations  at  similar  strain  rates  and  concentrations. 

On  the  basis  of  the  flows  so  far  Investigated,  it  is  possible  to  speculate  on  what 
would  happen  if  concentration  and  strain  rate  were  both  substantially  Increased.  It  is 
believed  that  the  flow  would  remain  suppressed  but  the  level  of  u''*'  would  Increase.  If 
this  should  be  the  case,  then  perhaps  the  very  substantial  difference  between  the  u'*^ 
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maxliui  measured  by  Relschman  and  Tlederman  (ref.  8)  and  those  measured  by  Rudd  (ref.  5) 
could  be  explained. 

Of  the  other  Investigations  which  obtained  extensive  velocity  measurements 
(using  LDA)  in  the  near  wall  region  of  a turbulent  drag  reducing  flow,  Relschman  and 
Tlederman  (ref.  8)  operated  at  essentially  the  same  range  of  strain  rates  as  we  did. 
Rumor  and  Sylvestor  (ref.  7)  were  very  slightly  higher,  and  Rudd  (ref.  5)  was  nearly  an 
order  of  magnitude  higher  (~10^  sec'^) . (Logan  [ref.  bj  did  not  report  strain  or  u • ) 

All  of  these  Investigators  used  100  ppm  or  more  of  a similar  or  longer  more  effective 
drag  reducing  polymer  than  WSR-301.  Thus,  we  expect  that  each  of  these  investigations 
measured  flows  for  which  t^e  near  wall  layers  could  be  characterized  as  suppressed 
turbulence  and  that  the  u*  measurements  reported  were  composed  mainly  of  polymer 
Induced  fluctuations  and  fluctuations  Imposed  by  the  outer  flow  which  were  not  Involved 
in  mixing  near  the  wall.  Since  the  water  table  experiments  reported  here  extended  only 
to  y of  100,  direct  comparison  Is  difficult  and  somewhat  conjectural.  Closed 
comparisons  would  have  been  possible  if  velocity  spectral  data  had  been  measured  by  the 
other  investigators. 

The  selective  seeding  of  the  flow  with  hydrogen  bubbles,  which  allowed  detailed 
spatial  resolution  and  spectral  measurement  of  the  velocity  to  be  made.  Is  believed  to 
be  novel.  Its  use  In  conjunction  with  visual  observations  of  the  free  surface 
disturbances  enabled  the  lamlnar-llke  polymer  Induced  velocity  fluctuation  flow  to  be 
identified. 

Comparisons  with  other  research  have  been  made  on  the  basis  of  wall  strain  rates. 

If  one  were  to  determine  Reynolds  numbers  appropriate  for  water  table  flows,  using 
thickness  as  the  length  parameter,  they  would  certainly  be  much  lower  than  the  Reynolds 
numbers  of  previous  investigations  (Refs.  5-9),  yet  the  very  suitability  of  drawing 
conclusions  from  a Reynolds  number  comparison  is  questionable  due  to  the  inherent 
differences  of  the  water  table  flow  from  those  in  pipes  or  over  a flat  plate.  Since  it 
is  fairly  well  established  that  drag  reduction  due  to  polymer  additives  is  a wall 
phenomenon,  using  wall  variables  appears  to  be  the  most  reasonable  approach  for 
comparison. 
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Fig.  2 Side  view  of  water  table  showing  Che  orientation  of  the  laser-Doppler 

anemometer.  The  detailed  view  depicts  the  selective  seeding  method  In  which 
bubbles  generated  by  electrolysis  and  introduced  at  a narrow  range  of  y are 
used  with  Che  LDA. 
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Fig.  3 Mean  velocity  of  water  flows,  U'  vs  y' 

laminar:  T > 3.04  ca/sec,  v ■ 3.25  cm/sec; 
turbulent : O u-^  > 3.42  cm/sec,  • u-^  ~ 3.68  cm/sec. 

Vertical  lines  Indicate  the  location  of  the  free  surface.  The  dashed  curves 
are  those  which  have  been  used  to  approximate  previously  published  data  for 
turbulent  boundary  layers  In  Newtonian  fluids. 
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laminor  noise  level 


Turbulent  water  u'  vs  y . All  u'  data  has  been  corrected  for  the  laminar 
LDA  noise  whose  uncorrected  level  Is  shown. 

O Ut  “ 3.42  cm/sec.  V ux  • 3.83  cm/sec. 


5 Water  turbulence  spectra;  as  for  all  spectral  data  presented  (Figs.  9,  12-15) 
full  scale  is  500  Hz  and  the  viirtlcal  amplitude  is  20  mV  rms  per  division, 
u^  ■ 3.71  cm/sec.  (a)  y • 58  (b)  y"*"  ■ 48  (c)  y " 14.  The  water  turbulence 
spectra  are  essentially  siitilar  to  these  over  the  range  of  u^  and  y+ 
Investigated. 
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Fig.  6 


Photographs  of  the  free  surface  when  turbulent  spots  have  been  triggered  in  an 
otherwise  undisturbed  flow  over  the  water  table.  The  heavy  lines  in  the  back- 
ground are  2.54  cm  apart  and  the  flow  depth  Is  typically  2 or  3 mm.  In  each 
case  the  fully  turbulent  depth  corresponds  to  u 3.6  cm/sec. 

(a)  Pure  water.  Over  the  range  studied  the  general  appearance,  size  and  growth 
rate  of  the  spot  are  essentially  Independent  of  depth  and  wall  shear  stress. 

(b)  11-13  ppm  WSR-301  solution.  The  wall  strain  rate  of  the  background 
laminar  flow  ( =1130  sec~^)  is  smaller  than  that  necessary  to  cause  polymer 
fluctuations.  The  turbulent  depth  is  2.8  mm  while  the  laminar  one  is  2.4  mm. 

(c)  11-13  ppm  WSR-301  solution.  The  wall  strain  rate  of  the  background  flow 
( =1210  aer“l)  high  enough  to  sustain  polymer  fluctuations  which  are 
visible  as  ripples  on  the  free  surface.  While  a spot  is  still  identifiable, 
it  is  much  less  distinct  than  in  (b).  The  turbulent  depth  is  2.0  mm  while  the 
laminar  one  is  1.9  mm. 
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Fig.  7 Observations  of  polymer  solution  flows  without  disturbance  as  a function  of 
concentration  and  strain  rate  (dU/dy).  O smooth  free  surface:  V rippled 
free  surface  indicating  the  presence  of  polymer  fluctuations  (cf.  the  back- 
ground of  Fig.  6c).  Fluctuations  are  expected  in  the  shaded  region. 


Fig.  8 Polymer  Induced  fluctuations,  u'  vs  y (in  rough  order  of  increasing  activity). 

(1)  lA  ppm,  u^  ••  3.71  cm/sec.  (2)  16  ppm,  u^  » 3.36  cm/sec. 

(3)  15  ppm,  u^  “ 3.65  cm/sec.  (4)  25  ppm,  u.^  ■ 3.73  cm/sec. 

(5)  32  ppm,  u^  “ 3.81  cm/sec.  The  dashed  curve  is  u'"*"  vs  y"*"  for  turbulent 

water.  A vertical  line  indicates  the  location  of  the  free  surface.  Dotted 
lines  indicate  extrapolation  beyond  data  points. 
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Fig.  9 Polymer  fluctuation  spectra:  16  ppm,  u^  « 3.36  cm/sec. 

(a)  y+  - 57  (b)  - 39  (c)  y+  - 12  (d)  y"^  - 8 

The  enriched  spectral  amplitude  at  50-150  Hz  is  clearly  visible  near  the  wall 
This  Is  the  same  flow  as  curve  2 in  Fig.  8 
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10  Polymer  Induced  fluctuation  flows,  U vs  y , roughly  In  order  of  Increasing 
polymer  activity. 

0 14  ppm,  ttj  ■ 3.71  cm/sec.  ^ 21  ppm,  u^  “ 3.72  cm/ sec. 

D IS  ppm,  > 3.65  cm/sec.  ^ 25  ppm,  u.^  - 3.73  cm/sec. 

O 16  ppm,  - 3.36  cm/sec.  ♦ 32  ppm,  • 3.85  cm/sec. 

These  lest  two  cases  are  at  the  asymptotic  limit  and  the  mean  velocity 
profiles  are  practically  identical  to  those  of  suppressed  turbulence  at  the 
sasie  u_  and  concentration. 


Fig.  11  Polymer  turbulence  u*  vs  y 

(1)  9 ppm,  u.^  ■ 3.72  cm/sec.  (2)  16  ppm,  i 
(3)  15  ppm,  • 3.25  cm/sec.  (4)  16  ppm,  i 
(5)  25  ppm,  * 3.73  cm/sec.  (6)  32  ppm,  i 
A vertical  line  indicates  the  location  of  the 
are  extrapolations  beyond  data  points. 


. - 3.52  cm/sec. 

> 3.85  cm/sec. 

. ■ 3.85  cm/sec. 

ree  surface,  and  dotted  lines 


Fig.  12  Polymer  turbulence  spectra;  9 ppm,  u.^  “ 3.72  cm/sec. 

(a)  y"*"  • 58  (b)  y"^  • 48  (c)  y'^  “ 14.  These  conditions  are  the  same  as 

presented  for  water  in  Fig.  5. 
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Fig.  13  Polymer  turbulence  spectra:  16  ppm,  u_  ■ 3.52  cm/sec. 

(a)  y'*'  - 41  (b)  y^  ■ 32  (c)  y ■ 13  (d)  y"^  ■ 8 ^ 

This  is  the  same  flow  as  for  curve  2 of  Fig.  11.  Note  that  at  y “ 
spectrum  is  similar  to  those  for  polymer  fluctuations.  (Fig.  9). 


13,  the 


Fig.  14 


Polymer  turbulence  spectra:  16  ppm,  u « 3.85  cm/sec. 
(a)  y+  - 54  (b)  y"*"  - 35  (c)  ‘ 


This  Is  the  same  flow  as  for  curve 


19  '^(d)  y"^  - 14. 
4 of  Fig.  11. 


Fig.  15  Polymer  turbulence  spectra:  25  ppm,  u -’3.73  cm/sec. 

(suppressed  turbulent  regime)  (a)  y ■ 50  (b)  y - 26  (c)  y - 12.  Tl 
Is  the  same  flow  as  curve  5 of  Fig.  11.  The  laminar  counterpart  of  this 
flow  had  essentially  identical  spectra. 
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Fig.  16  Mean  velocity  of  polymer  turbulence,  U va  y 
^ 9 ppm,  u ■ 3.72  cm/sec. 

O 16  ppm,  u^  * 3.52  cm/sec. 

• 16  ppm,  ■ 3.85  cm/sec. 

V lA  ppm,  u~  ■ 3.77  co/sec. 

0 21  ppm,  u]!  - 3.73  cm/sec. 

This  last  case  is  an  example  of  the  suppressed  turbulence  regime  and  is 
practically  identical  to  the  polymer  fluctuation  mean  velocity  profile 
under  the  same  conditions. 
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ABSTRACT 

A novel  method  of  improving  the  spatial 
resolution  of  a laser-Doppler  anemometer  in  water 
flows  has  been  developed.  Bubbles  generated  by 
electrolysis  from  a small  diameter  platinum  alloy  wire 
located  upstream  of  the  laser  scattering  voliae  is 
used  to  scatter  light  for  the  anemometer.  Most  of  the 
indigenous  scattering  particles  are  filtered  from  the 
flow.  The  bubble  generation  can  be  large  enough  to 
obtain  a nearly  continuous  velocity  signal  from  the 
laser  aoemometcr,  allowing  spectra  and  other  detailed 
velocity  measurements  to  be  obtained.  Results  are 
presented  comparing  hot-film  measurements  with  those 
obtained  using  the  selective  seeding  method. 
Comparisons  are  presented  for  a turbulent  flow  in  a 
fliae  and  for  laminar  and  turbulent  boundary  layer 
flows  involving  very  steep  velocity  gradients. 

NOMENCLATURE 

Cjj  drag  coefficient  of  circular  cylinder 

d diameter  of  bubble  wire  2Su 

B(22)  ratio  voltage  response  of  the  LDA  using  2^ 
beam  separation  to  that  of  the  hot  film  probe 
R(50)  ratio  voltage  raaponae  of  the  LEA  uaing  50am 
beam  aaparatlon  to  that  of  the  hot  film  probe 

“od 

Be  ■ Reynolds  maber  baaed  on  free  atraam 
valoclty  i wlra  diamatar 
U local  maaa  velocity  on  the  water  table 

U.  free  stream  valoclty 

P 

w local  bubbla  valoclty 

wall  friction  velocity 

a diataaea  from  bubbla  wire  to  acatteriag 

volume  in  the  stream  diraction 
y diataaea  meraal  to  the  wire  wake  or  distance 

away  from  tha  wall 

y*^  ■ ^ aon-dimanoioaal  diaunca  from  tha  wall 

V 

v blwatlc  viaeesity 


INTRODUCTION 

In  order  to  investigate  the  effects  of  strain 
rates  on  dilute  polymer  solutions,  it  was  necessary 
to  obtain  velocity  measurements  within  boundary 
layers  typically  thick  (Bertschy  t Abernathy,  1977). 
A laser-Doppler  anaiaoeter  (LDA)  was  chosen  to  obtain 
true  velocity  measurements  of  these  non-Newtonian 
flows;  however,  spatial  resolution  remained  a problem. 
The  scattering  voltae  formed  by  the  intersecting  laser 
beams  of  an  LDA  operated  in  the  dual  beam  mode  depends 
on  the  specific  geometry  of  the  optical  systmt  and  is 
typically  on  the  order  of  O.llmm  dia.  and  l.bn  in 
length.  If  the  long  axis  of  the  scattering  volume  was 
aligned  spanwise  to  the  flow,  aome  spatial  rasolution 
was  poasible  but  a sophisticated  optical  support  would 
have  been  necessary  to  traverse  the  LDA  system  to 
obtain  a flow  profile.  If  the  long  axis  was  aligned 
normally  to  the  surface,  across  the  velocity  gradient, 
then  apatlal  resolution  was  impossible.  Consequently 
for  shear  layers  of  the  dimensions  described,  it  was 
difficult  if  net  impossible  to  distinguish  the  signal 
of  a laminar  flow  from  one  of  turbulent  flow.  To 
ovorcome  this  limitation  in  spatial  resolution  of  an 
LDA  system,  a salactlve  seeding  technique  was 
davalopad  which  works  wall  with  water  flews. 

The  assance  of  the  technique  la  to  filter  out 
most  of  the  scattering  partlclas  from  the  water,  and 
than  introduce  a thin  shaat  of  hydrogen  or  oxygen 
bubbles  upstream  of  the  acatteriag  voliae,  and 
parpandlcular  to  its  long  axis.  Tha  latoraactlon  of 
tha  shaat  of  bubblaa  and  the  LDA  optical  veliaa  titan 
comatltutas  the  affoctlva  acattaring  volimM,  tkaraby 
raduciag  its  affactlva  laagth.  With  eufficiant  bubble 
density,  a maarly  continuous  valoclty  manaurMaat  can 
ba  obtaimad  which  is  aosantlal  for  noma  maaauraMnts. 

Tha  anphasia  of  this  paper  is  to  aatablish  the 
baala  for  uaiaR  tha  oalaetiva  naadlng  tachaiqua  and 
earring  its  raaults  with  those  abtalmad  from  standard 
laatriMntntien.  To  this  and,  maasuromaats  takaa  in  an 


open  channel  flow  with  an  LDA  eye ten  ualng  etandard 
natural  aeadlng  and  with  a hot-flln  aensor  were 
compared  to  ■aaauraments  obtained  using  selective 
seeding.  To  establish  the  validity  of  the  selective 
seeding  technique  in  a flow  with  steep  velocity 
gradients,  its  results  were  compared  with  those 
obtained  from  a hot  film  sensor  in  a flow  over  a 
water  uble. 


BUBBLE  VELOCITY 

Hydrogen  bubbles  generated  by  electrolysis  in 
water  have  been  used  for  generating  fluid  flow  markers; 
most  of  the  films  produced  under  the  direction  of  the 
National  Connittee  for  Fluid  Mechanics  Films  (Shapiro 
(1972))  used  the  bubble  technique  for  flow  visualisa- 
tion. Schraub  at.  al.  (1965)  were  early  developers  of 
the  technique  and  reported  data  suggesting  that  the 
generating  ifire  wake  had  negligible  effects  on  bubble 
velocity  at  a distance,  x,  greater  than  70  wire 
diameters.  More  recent  measurements  by  Grass  (1971) 
suggest  that  even  at  750  wire  diameters  the  wire  %iake 
causes  the  bubbles  to  move  measurably  more  slowly  than 
the  flow  outside  of  the  wake.  Since  each  of  these 
observations  was  derived  from  high  speed  motion 
pictures  and  consequently  had  substantial  uncertainty 
in  velocity,  we  undertook  an  independent  measurement 
of  bubble  velocity  in  the  generating  wire  wake.  The 
measurements  were  made  in  free  surface  water  flimM, 

60  channel  widths  from  the  entrance  where  flow  condi- 
tions arc  independent  of  distance  along  the  channel, 
er  in  the  discharge  of  a 2.54cm  nessla.  The  reaults 
are  shown  in  figure  1,  along  with  those  of  Crass  (1971) 
and  the  calculated  canterline  velocity  of  a laminar 
wake.  The  velocity  of  the  bubbles  u,  and  the  free 
stream  velocity  ware  both  obtaimad  from  LM 
maasuraments.  The  free  stream  velocity  mms  maaaurad 
without  a bubble  wire  la  place,  mslag  the  ambient 
particles  In  the  water  as  light  scattarars.  The 
results  are  for  four  dlifaraat  Baymolds  aMbars 
spaanlag  the  raage  of  free  stream  valecltlaa  genarally 
anceuatared  la  free  smrfaco  flows.  la  all  caaas,  the 
wire  dlaaater,  d,  mms  liv.  There  dees  met  appear  te 
be  any  Beymelda  wbsr  dap amd aura  far  the  valaclty 
defect,  idilcb  at  flrat  aeamad  eurprlalng.  Tha  dashed 
lias  dallad  the  lamlaar  wshe  la  a raaeenabla  fit  te 
tha  bubble  voleclty  data  far  900  < m/d  < 1000,  a 
eaawealeat  dlstaaea  raaga  far  the  ealeetlva  eaedlag 
tmehmlqma.  The  1 ami  mar  mmha  curve  is  the  theeretlcal 
eentar  line  vmloeitr  dmfaet  far  the  laminar  umke 
hahlad  a cue  dimaasienal  ehjoct  with  me  camaiidarmtiem 
of  bubble  gamaracien.  The  tbaeretleal  fen  ef  omeb  a 


FIGURE  1:  BUBBLE  VELOCITY  DEFECT 


UNCERTAINTY  SHOWN  IS  ASSOCIATED  SOLELY  WITH  ELECTRONICS 
O Re  • 39  T Re  - 7.1 

• Re  • 12  ^ Re  • 2.8 

LAMINAR  WAKE  CENTERLINE  DEFECT  (EQU.  3) 

GRASS  (1971) 


"o  4/r  V) 

and  la  given  la  many  rafaraaces  (e.g.  Schllchtlng, 
1968).  Tha  Rayaolds  mmber  dependence  appears 
eapllcltly  la  the  square  root  tarn  and  Implicitly  In 
tha  drag  ceafflclent,  C^.  A plot  of  versus 

6i^ 

Raynolds  nnbar  Is  shown  la  figure  2 where  the  values 
ef  Cji  as  a function  ef  ta  are  from  Trltton's  (1959) 
mnrrlaenul  maasurnaats  ualng  cylindrical  reds. 


noou  2.  m MMUCAL  ootmciBiT  or  m umhur 
auB  BoumoN  (vjo.  i)  roi  a cncoLa  ctuider, 
CptSa/A/r,  va  M.  Cg  «U  DBmMins  OBIMG 

WTTON'I  (19S9)  BATA. 


When  3 < Re  < 40,  the  product  only  varies  from 

4/ir 

1.26  to  1.3S.  We  have  taken  this  product  to  be  a 


constant  1.30  yielding  the  following  equation  for  the 
centerline  velocity  defect 


This  is  shown  as  a dashed  curve  in  figure  1 which  is  a 
reasonable  fit  to  the  measured  data  for  2.8  < Re  < 39. 
The  solid  curve  of  figure  1, 


-0.729 


due  to  Grass  (1971)  represents  a best  fit  to  his 
motion  picture  data.  The  uncertainty  in  his  measure- 
ments of  the  flow  Reynolds  number  was  not  given.  The 
uncertainty  displayed  in  figure  1 is  merely  associated 
with  the  difference  between  two  velocities  of  nearly 
equal  size,  each  known  to  ^ 11.  Greater  effort  was 
not  spent  in  reducing  the  uncertainty  since  the 
correction  Itself  is  small,  for  example 
6.5  t 2X  at  x/d  - 400. 

The  conclusion  to  be  drawn  from  figures  1 and  2 
is  that  hydrogen  bubbles  generated  by  a wire  in  the 
flow  do  follow  the  fluid  flow,  bearing  in  mind  that 
correction  is  necessary  for  the  wire  wake.  The 
correction  is  adequately  described  by  cqu.  2 for 
the  range  3 < Re  < 40. 


VEUtClTY  MEASUREMENTS 

The  velocity  measurements  of  figure  1 were 
presented  without  much  information  about  the  velocity 
signals  themselves.  Tyo  somewhat  different  techniques 
were  used.  Without  the  introduction  of  hydrogen 
bubbles,  one  must  rely  on  the  uniformly  distributed 
solid  psrticlss  in  the  wetsr  to  scsttsr  the  Iseer 
light.  Even  in  so-csllsd  filtered  tep  water  some 
solid  psrticlss  remain  to  scatter  light,  though  the 
particle  density  is  lew.  Using  a TSI  model  1090 
tracker  end  a 20mw  He-Me  laser  over  en  order  of 
nagnltude  variation  In  count  rate  was  possible  with 
the  roughly  filtered  meter  used,  depending  en  the 
gain  setting  of  the  phetedeteetor.  Typical  voltage 
signals  from  the  TSI  1090  tracker  ere  shown  in 
figure  3.  Each  oaclUoscepe  trace  represents 
1 second  of  dsu  obtained  from  the  discharge  of  the 
mettle.  The  time  average  velocity  met  nearly 
iadapeadaat  of  spetlsl  gradient  along  the  vertical 
nidplnns  end  equal  to  26ca/eoc.  Trace  (e)  is  low 
gala  200  eeuats/sec  (realitations  of  an  internally 


verified  velocity)  signal;  (b)  is  for  a high  gain 
setting  7000  counts/sec  signal;  and  (c)  is  for  the 
same  gain  setting  at  (a)  except  that  bubble  generating 
wire  was  located  400  wire  diameters  upstream  from  the 
LOA  scattering  volume.  The  hydrogen  bubbles  act  as 
very  effective  light  tcatterers  et  teen  by  compering 
traces  (a)  and  (c).  The  noise  of  the  signal  using 
selective  seeding,  (c),  it  lets  than  that  using  natural 
seeding  uniformly  distributed  simply  because  the 
selective  seeding  from  the  25i/  die.  wire  has  reduced 
the  effective  noasuring  volMe  size.  For  purposes  of 
averaging  and  spectral  analysis,  the  signal  obtained  by 


FIGURE  3:  OSCILLOSCOPE  TRACES  OF  VELOCITY  IN  THE 
FLUME  USING  LDA  AT  ABOUT  26  CM/SEC  WITH  (a)  NATURAL 
SEEDING.  200  VELOCITY  SAMPLES/SEC  WERE  MEASURED. 

(b)  NATURAL  SEEDING.  7000  VELOCITY  SAMPLES/SEC  WERE 
MEASURED  (c)  SELECTIVE  SEEDING.  THE  BUBBLE  WIRE  WAS 
AT  x/d  • 400  AMD  7000  VELOCITY  SAMPLES/SEC.  WERE 
MEASURED,  INDICATING  A VELOCITY  DEFECT  OF  ABOUT  lOX. 
WHEN  THE  BUBBLE  WIRE  WAS  TURNED  OFF,  200  SAMPLES/SEC. 
WERE  MEASURED.  ALL  TRACES  ARE  1 SEC,  DURATION. 


selective  sanding  can  be  treated  in  the  sane  manner 
as  the  signal  from  natural  seeding.  Since  the  count 
rate  is  higher,  it  is  often  possible  to  exanlns  higher 
frequency  fluctuations  in  velocity  using  selective 
seeding  than  when  using  the  natural  beckground 
pertleles  provided  the  hydrogen  bubble  can  actually 
follow  the  velocity  fluctuations  in  the  flow. 

VELOCITY  FLUCTUAnONS 

In  order  to  detemins  hew  well  the  hydrogen 
babbles  netunlly  follemod  the  vnleelcy  fluctuations  in 
the  flew,  the  spectral  mspenae  of  the  LBA  velocity 
signal  one  eonpared  with  that  obtalnad  fron  a eonven- 
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SIDE  VIEW  OF  MATER  TABLE  SHOWING  THE  ORIENTATION  OF  THE  LASER 


DOPPLER  ANEMOMETER  AND  DEPICTING  THE  SELECTIVE  SEEDING  METHOD 


FIGURE  4:  SCHDIAnC  OF  HATER  TABU  SBCWINC  THE  LDA  ARRANGEMENTS  AND  DEPICnHG  THE  8EUCTIVE  SEEDING  METHOD 


tieaal  bet-flla  enwoeeter.  TWe  different  flew 
eltuetlene  were  Inveetlfeted,  neaely,  thet  in  the 
15cn  wide,  free  eurfeee  eellbretlon  flime  end  thet  on 
e free  eurfeee  weter  teble  (flcure  A)  idiere  hlfh  etreln 
rete  flows  cen  be  eehleved  (the  velocity  typleelly 
verled  free  0 to  100  ea/eec  eeroee  the  eppraelaecelv 
Bh  depth).  In  the  eellbretlon  flime  the  nenlaim 
freqeeeey  of  velocity  flncteetlene  ie  on  the  order  of 
lOOb,  lAlle  en  the  enter  teble  fregweeclee  on  the 
order  of  MOHs  ere  obeerved.  It  eao  for  nee  on  the 
enter  teble  thet  the  eelectlve  rTtiflm  technlgne 
ene  first  developed. 


A conventions!  hot>fll*  snseoeeter  ees  need  ns  the 
stenderd  for  nessurlnc  velocity  fluctustlons.  The 
probe  wee  e TSI  nodel  1218  boundery  leyer  probe  (Slv 
dleaeter  end  qusrtt  costed  for  use  In  enter),  end 
ceoaseted  to  e TSI  nodel  1053B  constent  tsnpereture 
ensnoneter  bridge.  The  aon>llaser  reeponee  ees  then 
need  for  neesurenent.  The  constent  tsnpereture  eyeten 
hnd  e frequency  reeponee  well  into  the  kHt  rente,  well 
heynnd  the  freqetney  ef  ebeerved  flow  flMctustione.  In 


velocity  neesurint  eye tens,  the  voltete  entput  free 
the  LBA  treeher  or  tke  hot-filn  snplifler  ene  preeeoeed 


through  a Taktronlx  7LS  SpacCruB  Analyzer.  Such  a 
apactriBB  analyzer  aweepa  an  adjuatable  narrow  band 
digital  filter  acroas  a given  frequency  range.  The 
rate  of  awaep  of  the  filter  acroaa  the  range  is 
conatant.  hence  the  averaging  tlae  of  the  filtered 
algnal  near  zero  Hertz  la  leaa  than  the  period  of 
the  center  frequency.  A turbulent  velocity  algnal 
proceaaed  through  auch  an  Inatrvnent  will  naturally 
exhibit  conalderable  fluctuation  in  the  Intenaity 
at  the  zero  frequency  end  of  the  apectrun  which  ia  an 
artifact  of  the  electronlca  and  would  not  be  preaent 
if  the  aweep  rate  were  lower. 

The  7L5  Spectruo  Analyzer  waa  not  deaigned  to 
accept  a dc  algnal  froo  a high  Inpedance  aource  auch 
aa  the  TSl  LDA  ayateo.  Hence  it  waa  capacitor-coupled 
to  an  lapedance  Batching  operational  amplifier  of 
gain  1,  reaultlng  in  attenuation  of  the  algnal  below 
AHz.  Figure  S ahowa  the  output  of  the  Spectrum 
Analyzer  from  a dc  input  on  two  amplitude  acales 
uaed  in  the  data  preaented.  The  nolae  level  of  the 
ayataa  la  Iqw  and  the  output  ia  flat  abo\'e  lOHz.  The 
eacllloacope  output  automatically  dlaplaya  the  value 
of  four  adjuatable  paraaetera,  three  of  which  are 
uaeful  in  underatandiog  the  data  to  be  preaented.  At 


ncOU  5:  ttSFORSt  Of  7L5  SPIcnUM  ARALTm  TO  A DC 
(lOUL  AT  DirmnT  AMPUnCATIOH.  (a)  aoff/dlv. 
mriCAL  iriCTtAL  AKPUTDDE  (b)  S^/dir.  mnCAL 
iriCTIAL  AWLITUDE.  THE  BOEtZONTAL  SCAU  EONS  ntOM 
0 TO  SOOHa. 


the  top  center  la  the  vertical  intenaity  acale  in 
alllivolta  R.M.S.  per  centlaeter;  bottom  center  la 
the  width  of  the  adjuatable  digital  filter  in  Hertz; 
and  at  the  bottom  right  ia  the  abaclaaa  acale  in  Hertz 
per  centlaeter,  l.e.  per  grid  apace  in  the  photograph. 
All  data  preaented  waa  aaaaured  ualng  the  lOHz  filter 
and  covert  the  frequency  range  from  zero  to  500Hz. 

In  order  to  compare  aeanlngfully  the  LDA  output 
voltage  fluctuatlona  with  thoae  of  a TS1-1053B- 
eonatant  teaparature  hot-flla  anemometer,  it  it 
neceatary  to  calibrate  the  hot-flla  ayatam  to  have  a 
voltage  veraut  mean  velocity  curve  for  a given  over- 
heat ratio.  Thla  waa  done  in  the  conventional  way 
ualng  the  LDA  ayttem  with  natural  aeeding  aa  the 
velocity  atandard.  The  local  alope  of  the  voltage- 
velocity  curve  waa  uaad  at  the  linearized  reaponae  of 
the  hot-flla  ayatem.  The  LDA  ayatem  ia  linear  in  Ita 
reaponae;  however,  a ahift  in  aeparatlon  in  the 
optical  beam  aplltter  changet  the  elope  of  the 
voltage-velocity  curve.  By  aelectlng  a convenient 
mean  velocity  and  beam  aeparatlon  it  waa  poatlble  to 
adjuat  the  reaponae  ratio,  R(22)  or  R(50)  of  the  LDA 
output  voltage  to  the  hot-flla  output  voltage  to  be 


rieOEE  6;  moaTT  sfectra  in  the  fume  flow  at 

AO  CM/SBC.  (•)  IDA  (SOm  BEAM  SEFAEATIOH)  - 
HAIVSAL  IBIDIHC  (b)  KT-FZUI  tnSOE.  E(S0)  - A.  1 AND 
TD  OAIN  HAS  KIN  ADJDITD  OH  TK  7U  IT  A FACTOE  OF 
A TD  COWAU  HITH  te.  FOE  OOKAEISaN  HITE  FICUEE  7. 
E(22)  - 1,7. 


frcquancy. 

Spectra  of  the  laalnar  flow  velocity  al(nal  at 
■eaaured  by  the  aelactlve  teed log  LOA  aethod  and  by  a 
hot-flln  probe  are  coaparad  In  figure  9. 


riGURE  9:  VELOCITY  SEECTIIA  OF  A LAMINAR  UATEK  TABLE 
FLOW  (PICUKE  8).  (a)  LOA  (22),  y - O.ObbCM.y'*’  - 25 

(b)  LDA  (22),  y - 0.046CM,y'*‘  - 18.  WHEN  10  < y'*’  < 25 
FLOW  EXCITATION  OF  A NATURAL  FREQUENCY  OF  THE  BUBBU 
WIRE  SUPPORT  IS  OFTEN  OBSERVES  BY  HAVING  A PEAK  IN  THE 
SPECTRUM  AS  SHOWN  HERE  AT  ABOUT  120HZ. 

(e)  BOT-FIUt  SENSOR,  y - 0.051CH,y'*'  - 20. 

Ceaparlng  the  hot-flla  algnal  (flgura  9c)  with  either 
LOA  algnal  (figure  9a  or  9b)  ahowa  the  nolaa  level  of 
the  aelactlve  aeedlag  technique  In  auch  a flow. 

Laalnar  flow  apectra  at  y*  below  25  often  ashlblt  a 
eplka  in  Intenaity  in  the  range  of  100  to  200Ht.  Thla 
artifact  of  the  technique  will  be  dlacuaaed  In  the 
aeat  aectlon. 

Terbaleat  velocity  apectra  aaaaurad  by  the  two 
teehaigaea  are  proeaated  in  figure  10.  The  only 
slgalf leant  difference  between  the  hot-flln 


noaauranentt  and  the  aelectlvely  aeeded  LDA  neature- 
nente  la  the  larger  nolee  level  at  high  frequency. 

Proa  figure  9t  It  appears  that  the  LOA  nolte  la  broad- 
band and  In  fact  ahould  be  thought  of  at  underlying 
the  turbulent  tpeetT\n  (figure  lOt)  alto. 


FIGURE  10:  VELOCITY  SPECIRA  OF  A TURBULENT  HATER 
TABLE  FLOW  (FIGURE  8).  (a)  LDA  (22),  y » 0.147CM. 
y*  • bl.  (b)  HOT  FILM  SENSOR,  y • 0,155CM,  y'*'  • 67 
R(22)  - 1.9. 


Except  for  the  higher  broad-band  noiae  leval.  It 
hat  been  danonatrated  that  the  LDA  aelactlve  aeedlng 
technique  natchea  Che  perfotnance  of  the  hot-flla 
anew one try.  Bubbles  froe  the  2Su  die.  wire  do  follow 
the  fluctuatlona  in  velocity  and  neature  the  average 
velocity  whan  corrected  by  aqu.  (2)  for  the  wire  woke. 
Schraub  at.  al  (1965)  concluded  fren  calculatlona  chat 
bubblea  of  the  dlanatar  uaed  In  theae  flows  will  in 
fact  respond  to  fluctuations  in  the  kiloHerts  ra^e; 
hwever,  we  have  not  found  It  neceesary  to  experlaen- 
tally  verify  the  upper  Halt  of  the  reaponse  capability 
of  the  syataa. 

The  selective  aeedlng  technique  extends  the 
usefulness  of  the  LDA  to  flow  raglona  of  Urge  velocity 
gradients  and  even  reduces  the  nolaa  level  in  ceaea 
where  the  naan  gradient  is  aero  or  anall.  It  has  been 
lavelueble  to  the  authors  la  reaearch  lnvolvli«  dilute 
polynar  solutlena  lAera  the  eppUcablUty  of  the  hot- 


Analyser  for  the  LDA  was  selected  to  be  2 or  4 times 
that  of  the  hot-flln  ao  that  the  photographs  of  the 
spectra  are  nearly  directly  comparable. 

The  turbulent  velocity  spectrum  In  the  calibra- 
tion fliaie  at  40  cm/eec  velocity  is  shown  in  figure  6 
as  maasured  by  the  LOA  system  using  natural  aceding 
and  by  the  hot-film  anemometer.  The  comparison  of 
the  two  spectra  shows  the  noise  level  introduced 
above  150  Hz  by  the  LDA  system.  The  turbulent 
spectrum  in  the  fltsae  is  concentrated  primarily 
below  100  Hz. 

Figure  7 shows  a comparison  of  the  same  flow 
maasured  three  different  ways.  The  voltage 
sansltlvlcy  of  the  LOA  is  a factor  of  2.4  less  than 
it  was  in  figure  6 since  the  beam  separation  was 
reduced  from  SOmm  to  22mm  to  facilitate  later 


(a)  NATOBAL  SnDlW:  (b)  lOUU  WIU  AT  x/d  - 100 
(c)  BOIBU  WIRE  AT  s/d  - 300. 


A comparison  of  figure  7(a)  and  7(c)  shows  the  lower 
noise  level  of  the  selective  seeding  technique 
(trace  (c))  compared  to  the  natural  seeding  at  high 
frequencies.  In  addition,  trace  (c)  at  500  diameters 
behind  the  wire  more  nearly  mimics  the  very  low 
frequency  response  of  the  natural  seeding  (a)  without 
a wire  present  than  does  (b)  taken  only  100  diameters 
behind  the  bubble  generation  wire.  The  selective 
seeding  technique,  by  easily  allowing  a high  data 
count  rate  and  by  reducing  the  effective  size  of  the 
scattering  volume  reduces  the  high  frequency  noise 
level  and  approaches  the  response  of  the  hot-film 
probe.  (Compare  (6a  and  7c)). 

The  flow  in  the  flume  is  restricted  to  low 
frequencies;  the  flow  on  the  free  surface  water 
table  on  the  other  hand  has  higher  frequency 
components.  The  uncorrected  velocity  profiles  of  a 
laminar  and  a turbulent  flow  of  water  on  the  water 
table  are  shown  In  figure  8.  Both  of  these  profiles 


U lam/ssc) 

riCtniE  8:  LAMINAR  AND  TURBULENT  VELOCin  PROFILES  OF 
A SHEAR  FLOW  OVER  THE  HATER  TABLE  AS  MEASURED  BY  THE 
SELECTIVE  SEEDING  METHOD  WITH  s/d  - 400.  THE  SOLID 
SYMBOLS  ARE  THE  MEASURED  VELOCITY  OF  POWDER  SPRINKLED 
ON  THE  FREE  SURFACE.  THE  ARROW  AT  0.004CM  INDICATES 
THE  LOCATION  OF  THE  WAU  USING  EXTRAPOLATION  OF  THE 
SLOPE  DETERMINED  BY  THE  HALL  SHEAR  STRESS  AND  THE 
LOWER  DATA  POINTS. 


Comparison  with  standard  LOA  using  natural  aamdlng  is 
Impossible  alnce  the  velocity  gradient  is  too  eteep 
for  ouch  meesuramants.  The  bubble  wire  was  aligned 
epenrlee  to  the  flew  400  diameters  upstream  of  the 
optical  crossing  voliae.  Details  are  praaeatmd  in  the 
nest  section.  The  LDA  output  voltage  was  averaged  by 
a cootlnuously  Imtegratlag  op-amp.  with  an  accuracy  of 
O.U  from  de  to  the  ooveral  blloHartt  raiwe  In 


fils  ancoooctry  la  doubtful  at  bast. 

COMMENTS  ON  USING  THE  TECHNIQUE 

Expcrianca  Indlcatad  that  great  care  oust  be 
exercised  to  ■alntaln  good  quality  ■eaturementt  from 
the  electrolysis  bubbles.  So,  the  wire  was 
occasionally  claanad  using  a soft  artist's  paint 
brush  to  remove  accumulated  dirt  (80%  Ft  20%  Ir 
proved  more  durable  than  100%  Ft  and  was  supplied 
by  Sigmund  Cohn  Corp.,  Mt.  Vernon,  N.  Y.).  Between 
brushlngs  the  wire  was  kept  clean  by  changing  the 
polarity  frequently  and  slewing  the  voltage 
(80v  to  300v)  In  order  to  drive  off  large  gas  bubbles. 
(As  a side  note,  there  were  also  unusual  difficulties 
in  using  the  hot-film  sensor  since  the  water  could 
not  be  de-lonlzed.  Ion  rich  water  was  needed  to 
t generate  good  electrolysis  bubbles.  To  prevent  severe 

* drift,  the  sensor  had  to  be  brushed  before  each 

raading.) 

When  using  the  selective  seeding  technique  In  the 
environment  of  a high  shear  flow  there  are  several 
I points  to  kaep  In  mind.  Tha  mean  velocity  profiles  of 

figure  8 were  obtained  by  moving  the  LDA  system  only 
i once.  Instead  the  bubble  wire  was  movad,  alimlnatlng 

the  oaed  for  a sophisticated  LDA  mount  or  the  need  to 
make  continual  optical  adjustments.  To  facilitate 
data  taking  naar  a boundary,  the  optical  scattering 
voltae  must  be  partially  burlad  in  the  glass  wall. 

Such  meaauraments  with  selective  seeding  need  no 
significant  corraction  since  the  effective  scattering 
voliac  la  always  completely  within  the  flow.  The  data 
of  figure  8 taken  naarest  to  the  wall  show  a definite 
tendency  toward  spuriously  high  velocities.  In  all 
likelihood,  this  was  caused  by  driving  the  bubble  wire 
into  the  wall.  Since  tha  wire  support  was  flexible, 
the  wire  rose  into  the  flew  as  the  support  was  driven 
Into  the  wall.  Hhila  relative  wire  position  was  very 
accurately  known  (using  a dial  indicator  on  the 
euppert),  the  absolute  elevation  above  the  wall  could 
only  be  roughly  determined  by  finding  the  location  of 
minimum  velocity  raading.  Once  e profile  was  Obtained 
a better  determination  of  the  wall  location  was 
eccompllahed  by  extrapolating  the  profile  to  aero 
velocity,  dlacounting  thoae  polnta  which  aeamed 
apurieualy  high.  For  the  water  table  data  the  fluid 
depth  gave  an  independent  nethod  of  detemialng  the 
wall  ahaer  etraaa  aad  hence  the  velocity  profile* a 
elope  at  the  wall,  tinea  the  accuracy  of  thia 
■aeawrMent  waa  high,  it  wma  need  to  aid  in  locating 
the  wnll.  la  figure  I the  poaition  of  the  nmll  wma 


found  to  be  at  y : 0.004  cm  and  a small  arrow  Indicates 
the  location.  Such  an  Independent  knowledge  of  the 
wall  shear  stress  Is  frequently  obtainable. 

Another  phenomenon  related  to  the  wire  support 
can  be  seen  In  the  spectrum  of  figure  9b.  When  ‘ 

maasurlng  In  laminar  flows  near  the  wall,  discrete 
frequencies  were  often  observed  in  the  velocity 
spectrum.  These  appear  to  be  caused  by  flow 
excitation  of  a natural  frequency  of  the  wire  support. 

For  the  laminar  profile  of  figure  8 such  a spike  was 
observed  when  10  < y < 25. 

Further  Investigation  of  figure  8 shows  that  as 
the  wire  nears  the  free  surface,  the  velocity  measure- 
ment Is  spuriously  low.  This  is  probably  the  result 
of  Che  wire  incerectlng  with  the  surface  primarily  due 
to  surface  tension  effects. 

To  obtain  a surface  velocity  independently  of 

/ 

selective  seeding  with  electrolysis  bubbles,  talcum 
po%<der  was  sprinkled  on  the  free  surface  and  used  as 
a prefarentlal  scattaring  agent  with  the  LDA.  While 
high  density  readings  were  not  possible,  approximate 
mean  valocities  were  obtainable  and  are  depicted  es 
solid  data  points  In  figure  8.  They  appear  to  run 
about  11%  higher  than  the  extrapolation  from  the 
selective  seeding  data.  From  the  analysis  of  the 
bubble  velocity  defect,  a difference  of  ebout  6.5% 
from  true  free  stream  velocity  wes  predicted.  Since 
talciB  powder  la  denser  than  water  (sp.  gr.  2. 6-2. 8) 

It  will  accelerate  to  a higher  velocity  at  it  moves 
down  the  Incline  of  the  water  table  and  through  the 
interaction  of  the  film  of  talcisD  powder  which 

forms  on  the  surface,  free  surface  boundary  layers  I 

were  expected.  When  lycopodiut  poivder  | 

(density^  water)  was  used,  the  surface  velocity  i 

meaaurament  ran  1 or  2%  less  than  when  using  talcum  j 

powder.  j 
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